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We examine the motion and tidal dynamics of a nonrotating black hole placed within a post- 
Newtonian external spacetime. The black hole's gravity is described accurately to all orders in 
Gm/c 2 r, where m is the black-hole mass and r is the distance to the black hole. The tidal pertur- 
bation created by the external environment is treated as a small perturbation. At a large distance 
from the black hole, the gravitational field of the external distribution of matter is assumed to be 
sufficiently weak to be adequately described by the (first) post-Newtonian approximation to general 
relativity. There, the black hole is treated as a monopole contribution to the total gravitational 
field. There exists an overlap in the domains of validity of each description, and the black-hole 
and post-Newtonian metrics are matched in the overlap. The matching procedure produces (i) a 
justification of the statement that a nonrotating black hole is a post-Newtonian monopole; (ii) a 
complete characterization of the coordinate transformation between the inertial, barycentric frame 
and the accelerated, black-hole frame; (iii) the equations of motion for the black hole; and (iv) the 
gravito-electric and gravito-magnetic tidal fields acting on the black hole. We first calculate the 
equations of motion and tidal fields by making no assumptions regarding the nature of the post- 
Newtonian environment; this could contain a continuous distribution of matter (so as to model a 
galactic core) or any number of condensed bodies. We next specialize our discussion to a situation 
in which the black hole is a member of a post-Newtonian two-body system. As an application of 
our results, we examine the geometry of the deformed event horizon and calculate the tidal heating 
of the black hole, the rate at which it acquires mass as a result of its tidal interaction with the 
companion body. 



I. INTRODUCTION AND SUMMARY 

A. This work and its context 

How does a black hole move in an external spacetime, 
and what effects do the tidal fields created in the external 
spacetime have on the black hole? These are the ques- 
tions that are investigated in this work, in a context in 
which the black hole is nonrotating and the gravity of the 
external universe is sufficiently weak to be adequately de- 
scribed by the post-Newtonian approximation to general 
relativity. This work is a continuation of a line of inquiry 
that was initiated by Manasse pQ in the early nineteen 
sixties, and that has been pursued to the present day. 

To pose our questions more precisely, and to better 
discuss the place of this work in the context of what was 
achieved previously, we introduce two length scales that 
are relevant to this problem. The first is set by m. the 
mass of the black hole, which gives rise to an associated 
length scale M := Gm/c 2 , the gravitational radius of 
the black hole. The second is 1Z, the radius of curvature 
of the external spacetime, evaluated at the black hole's 
position. Our work, and all others that preceded it, is 
carried out in a context in which M/7Z < 1, so that 
there is a clean separation between these scales. Only in 
this context can one meaningfully speak of a black hole 
moving in an external spacetime; when M is comparable 
to 7Z, no distinction can be made between the "black 
hole" and the "external spacetime." 

As a concrete example we may consider a situation in 
which the black hole is a member of a binary system. 
Then 7Z ~ \/b 3 /Mt t, where b is the separation between 



the bodies and M to t ■— G(m + m') /c 2 is a measure of the 
total mass within the system [m! is the external mass). 
In this case we have 

M M / Mtot \ a/" 
TZ ~ M tot V b J ' 

and for this work, this is required to be small. There 
are two particular ways to achieve this. In the small- 
hole approximation the black-hole mass is assumed to be 
much smaller than the external mass, so that M/M tot ~ 
m/m' <C 1; then M/1Z is small irrespective of the size of 
Mtot/b, and the binary system can be strongly relativis- 
tic. In the weak- field approximation it is M tot /b that is 
assumed to be small, while the mass ratio is left uncon- 
strained. 

Our work is concerned with the weak-field approxima- 
tion. The black hole is placed within a post-Newtonian 
external spacetime, and the external gravitational po- 
tentials determine its motion as well as the tidal gravity 
acting upon it. We determine the motion of the black 
hole, the tidal fields, and the effects of the tidal fields 
on the structure of spacetime around the black hole, all 
within the post-Newtonian approximation to general rel- 
ativity. At first, we do not specify the nature of the post- 
Newtonian environment. We leave it completely general; 
the black hole might be immersed within a smooth dis- 
tribution of matter (a model for a galactic core, for ex- 
ample), or it might be part of an A^-body system (with 
the number, nature, and state of motion of the bodies 
left arbitrary). As our work progresses, we specialize our 
results to a two-body system undergoing generic orbital 
motion, and finally we examine the special case of circu- 
lar orbits. 
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The motion of a black hole in an arbitrary external 
spacetime was first investigated by D'Eath j2j |3j 0] and 
Kates [5], who showed that in the limit M/1Z — > 0, the 
black hole moves on a geodesic of the external space- 
time. In this limit the black hole behaves as a test mass, 
in spite of the fact that the self-gravity of the black hole 
never ceases to be strong. The corrections to geodesic 
motion produced by the coupling of the black-hole spin 
with the curvature of the external spacetime were worked 
out by Thorne and Hartle [6] , who also obtained preces- 
sion equations for the spin vector. These authors ex- 
ploited the power of matched asymptotic expansions in 
their derivation of the equations of motion. In their ap- 
proach, the metric of the black hole (deformed by the 
conditions in the external spacetime) is matched to the 
metric of the external spacetime (perturbed by the mov- 
ing black hole) . The matching is carried out in a region in 
which both descriptions are valid, and it produces both 
the equations of motion and the tidal fields, with only the 
Einstein field equations as additional input. Our work is 
a continuation of this program. 

These investigations were next specialized to systems 
for which the gravity of the external spacetime is weak; 
this is the context that interests us in this paper. Dcmi- 
anski and Grishchuk [7] showed that to leading order in 
a post-Newtonian expansion of the external gravity, the 
black hole moves according to the Newtonian equations 
of motion. Their results were generalized to first post- 
Newtonian order by D'Eath [3] and Damour [8 , who 
found agreement between the equations of motion for 
black holes in binary systems and the standard (Einstein- 
Infcld-Hoffman) equations of motion of post-Newtonian 
theory. Our work is a continuation of this effort, and our 
results are slightly more general than theirs: While our 
black hole is still immersed within a post-Newtonian en- 
vironment, this environment is completely general, and 
the black hole is not required to be a member of a bi- 
nary system. When, however, we specialize our results 
to this particular case, we recover the results of D'Eath 
and Damour. 

The motion of a black hole in a post-Newtonian ex- 
ternal spacetime is well understood, and our contribu- 
tion to this understanding is a relatively minor one. The 
same cannot be said, however, of the effects of the ex- 
ternal tidal gravity on the black hole, which have not 
been much discussed in the literature. This is the true 
focus of this work, and our main goal in this paper is to 
calculate the post-Newtonian tidal fields acting on the 
black hole, and to explore the physical consequences of 
the tidal interaction. 

We are not claiming that ours is the first calculation 
of post-Newtonian tidal fields acting on a self-gravitating 
body. It is not. In their pioneering work on relativistic 
celestial mechanics, Damour, Soffel, and Xu [5J [TUJ QT] 
calculated the post-Newtonian tidal fields acting on 
an arbitrarily-structured body with weak internal grav- 
ity. This work was recently generalized to arbitrarily- 
structured, strongly-gravitating bodies by Racine and 



Flanagan [T^] . Our work is concerned instead with a very 
specific type of strongly self-gravitating body: a non- 
rotating black hole. We calculate the post-Newtonian 
tidal fields acting on this black hole, and observe that 
they are the same as those obtained by Damour, Soffel, 
and Xu in the case of weakly self-gravitating monopolcs. 
We confirm, therefore, the general expectation (known 
as the "effacement principle") that the post-Newtonian 
tidal fields must depend on the body's multipole mo- 
ments only (in addition to the conditions in the external 
spacetime), and not on additional details concerning its 
internal structure. 

The effects of tidal fields on the structure of space- 
time around a black hole were first investigated by Man- 
asse PP, who provided an essential input to the work 
later carried out by D'Eath, Kates, Thorne, and Har- 
tle. Adopting the small-hole approximation defined pre- 
viously, Manasse calculated the metric around a small 
black hole that falls radially toward a much larger black 
hole. Each black hole was taken to be nonrotating, and 
the small hole was taken to move on a geodesic of the 
(unperturbed) Schwarzschild spacetime of the large hole. 
The tidal gravity exerted by the large black hole cre- 
ates a perturbation in the Schwarzschild metric of the 
small hole, and employing the techniques of Regge and 
Wheeler |T3] , Manasse was able to calculate this pertur- 
bation in a local neighborhood of the black hole. His 
metric is expressed as an expansion in powers of r/lZ, 
where r is the distance to the small hole; it is accurate 
through second order in r/lZ, and it is valid to all or- 
ders in M/r, which measures the strength of the small 
hole's self-gravity. The case of circular motion around a 
large Schwarzschild black hole was treated much later by 
Poisson P3]. 

The methods used by Manasse are not necessarily re- 
stricted to the small- hole approximation. Alvi jT5j US] 
realized that these methods could be seamlessly ex- 
tended to the general setting defined by the requirement 
M/TZ <C 1, which includes both the small-hole and weak- 
field approximations as special cases. Alvi exploited this 
insight to calculate the tidal fields acting on a black hole 
in a post-Newtonian binary system (perhaps with an- 
other black hole). In Alvi's work, the two bodies have 
comparable masses and the black hole has a significant 
influence on the geometry of the external spacetime. Alvi 
calculated the tidal fields to leading (Newtonian) order in 
the post-Newtonian approximation to general relativity, 
for circular orbits. The metric of the distorted black hole 
was next joined to the post-Newtonian two-body metric, 
and the global metric was presented in a single coordi- 
nate system that corotates with the system. (In Alvi's 
original work there is a discontinuity in the metric at the 
common boundary between the two descriptions. The 
joint was made continuous in a follow-up paper by Yunes 
et al. QI].) 

Alvi's insight was exploited by Poisson [18j in a cal- 
culation of the metric of a nonrotating black hole placed 
within an arbitrary tidal environment (still restricted by 



M/1Z <C 1). Working, like Manasse [T], in a local neigh- 
borhood of the black hole, Poisson was able to calculate 
the metric through third order in r/TZ, while keeping the 
expressions accurate to all orders in M/r. Poisson's met- 
ric is parameterized by a number of tidal moments, freely- 
specifiable tensorial functions of time that characterize 
the black hole's tidal environment. This metric gives a 
general description of the spacetime around a black hole 
in any tidal environment, but a more complete descrip- 
tion requires the determination of the tidal moments. 
This is our task in this paper: We calculate the tidal 
moments for the specific case described above, in which 
the black hole is immersed within a post-Newtonian ex- 
ternal spacetime. This is a generalization of Alvi's work 
[ISl [T5J [T7] : We appeal to the weak-field approximation, 
we calculate the tidal fields created by an arbitrary post- 
Newtonian spacetime, and we do so to a higher order of 
accuracy than what was achieved by Alvi. 



B. Our results 

The metric of a nonrotating black hole immersed in 
a tidal environment is expressed as a perturbation of 
the Schwarzschild metric. We take the black hole to 
have a mass m when it is in complete isolation (un- 
perturbed), and we denote its gravitational radius by 
M := Gm/c 2 . The strength of the tidal perturbation 
is measured by the inverse length scale TZ^ 1 , and we as- 
sume that M/TZ <C 1; the tidal perturbation is weak. 
In addition, we assume that the black hole moves in an 
empty region of spacetime, so that in the hole's neighbor- 
hood SB, the perturbed metric satisfies the vacuum field 
equations linearized about the exact Schwarzschild solu- 
tion. We present the metric in the comoving reference 
frame of the black hole, in a quasi-Cartesian system of 
coordinates x a = (x°,x a ) = (ci,x,y,z) that enforce the 
harmonic conditions dp(y/—gg al3 ) = 0. Our convention 
is that Greek indices run from to 3, while Latin indices 
cover the spatial coordinates and run from 1 to 3. We 
raise and lower Latin indices with the Euclidean metric 
5 a b, and we let e a bc denote the permutation symbol of 
ordinary vector calculus (with £123 = e xyz = 1). 

The time-time and time-space components of the 
black-hole metric are (Sec. Ill) 



.9oo 



1 - M/f 1 
~ 1 + M/f ~ 'c 2 

+0(f 3 /1Z 3 ), 
2 

+0(f 3 /lZ 3 ) 



(1 - M/f) 2 £ ab (t)x a x b 



(1.1) 



(1 - M/f)(l + M/f) 2 e abp B p c (t)x b x c 



(1.2) 



for our purposes here we shall not need an expression 
for g a i, the space-space components of the metric. The 
metric is expressed as an expansion in powers of f/TZ, 
the ratio of f := \/x 2 + y 2 + z 2 , the distance from the 
black hole, to 1Z, the local radius of curvature of the 




FIG. 1: The post-Newtonian domain 3), the black-hole neigh- 
borhood SB, and the overlap region 6 '. The post-Newtonian 
domain is depicted in blue (light grey), and it includes the 
green (dark grey) annulus that surrounds the black hole. The 
black-hole neighborhood is drawn as a yellow (white) disk 
around the black hole, and it also includes the green (dark 
grey) annulus. The black hole is represented as a black disk. 
The overlap region is the union of SB and Q> that is shown in 
green (dark grey), the union of the blue and yellow colors. 



external spacetime. The metric is valid in the black hole's 



local neighborhood S3, which is defined by r < 



(see 



Fig. 1); we demand that r max /TZ be small, but within S3 
the ratio M/f is allowed to be arbitrarily large. 



The first term on the right-hand side of Eq. ( 1.1 1 is 



the Schwarzschild piece of the metric; it is expressed in 
harmonic coordinates, and in these coordinates the (un- 
perturbed) event horizon is situated at f = M. The sec- 



ond term and the right-hand side of Eq. (1.2) represent 



the tidal perturbation. The tensorial functions £ a b(i) and 
B a b(t) are the tidal moments (Sec. II), and it is these ten- 
sors that characterize the black hole's tidal environment. 
The tidal moments are symmetric and tracefree (STF) 
tensors, in the sense that Ef, a = £ a b and 8 ab £ a b = 0, with 
similar relations holding for B ab - The tensors c~ 2 £ a b and 
c~ 3 B a b have a dimension of squared inverse length, and 
their scale defines 1Z, the local radius of curvature; we 
have c~ 2 £ ab ~ 1Z~ 2 and c~ 3 B ab ~ 1Z~ 2 . The tidal mo- 
ments are not determined by solving the Einstein field 
equations in SB. They are a priori arbitrary, and their 
determination is accomplished by matching the black- 
hole metric to a global metric defined in a domain that 
is much larger than SB. In this work the tidal moments 
are determined by placing the black hole within a post- 
Newtonian environment; the global metric is obtained 
by solving the Einstein field equations in the first post- 
Newtonian approximation. 

We next describe this post-Newtonian environment. 
We consider a spatial domain £F that is much larger than 
SB, the black-hole neighborhood (see Fig. 1). This do- 
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main contains an arbitrary distribution of matter 1 , and 
it is assumed that everywhere within Qi, gravity is suf- 
ficiently weak to be adequately described by the post- 
Newtonian approximation to general relativity. The do- 
main is spatially limited by a sphere of radius r noar cen- 
tered on the post-Newtonian bary center. This sphere 
marks the boundary of the near zone: If T is a typi- 
cal time scale for processes taking place within S, and 
if A c = cT is a typical wavelength of the gravitational 
waves escaping the domain, then r noar < A c . The do- 
main also excludes a sphere of radius r lmn centered on 
the black hole, inside which gravity is too strong to be 
adequately described by post-Newtonian theory. We de- 
mand both that M/f min -C 1 and f min /TZ <C 1, which 
is possible when M <C 7Z. There exists an overlap re- 
gion G between the black-hole neighborhood 3$ and the 
post-Newtonian domain S . This region is described by 
Praia < T < f m ax, and we assume that there is no matter 
in G . So while matter is present somewhere within S>, we 
assume that the black hole is moving in an empty region 
of spacetimc. 

The metric in S> can be expressed as a post- Newtonian 
expansion of the form (Sec. IV) 



5oo 

90a 
9ab 



-l + ^U+^-U 2 ) + 0(c- 6 ), (1.3) 



rU a + 0(c- 5 ), 



l + ^U)S ab + 0(c- i ), 



(1.4) 
(1.5) 



in which U is a Newtonian potential, U a a vector po- 
tential, and W a post-Newtonian potential; the metric 
is presented in harmonic coordinates x a = (x°,x a ) = 
(ct, x, y, z). These barycentric coordinates differ from the 
black-hole coordinates (ci, x, y, z) introduced previously; 
the black-hole and post-Newtonian metrics are presented 
in different coordinate systems. In spite of the fact that 
each system is harmonic, the coordinates are indeed dis- 
tinct: They are defined in different domains (x Q in 3/, x a 
in 3$) , and they have a different spatial origin (x a is cen- 
tered on the post-Newtonian barycenter, whose position 
is fixed in the global reference frame, while i° is centered 
on the moving black hole). 

In the overlap region G the spacetime is empty of mat- 
ter, and the potentials U, U a , and satisfy the vacuum 
field equations of post-Newtonian theory. The Newto- 



nian potential, for example, must satisfy Laplace's equa- 
tion in flat space, V 2 C/ = 0. The solution must account 
for the presence of a black hole, and it must also account 
for the presence of matter outside G. We treat the black 
hole as a post-Newtonian monopole, and we write 



U(t,x) 



Gm 



\x-z(t)\ 



+ E/cxt 



(1.6) 



in which the three-dimensional vector z(t) denotes the 
position of the black hole in the barycentric coordi- 
nates. The external potential C/ ext is created by the mat- 
ter outside G, and within the overlap region we have 
V 2 £/ext = 0. The potentials U a and "J are handled in a 
similar fashion (Sec. IV), and in this way we construct 
the post-Newtonian metric in G. 



The post-Newtonian metric of Eqs. (1.3 1— (1.5 1 and the 



black-hole metric of Eqs. ( 1.1 )— ( 1.2 ) both give a valid 



description of the gravitational field in G . The metrics 
must agree in the overlap region, and matching them de- 
termines the equations of motion for z(t) as well as the 
tidal moments £ a b(i) and B ab (t). This matching, how- 
ever, can only be done after the post-Newtonian metric 
is transformed from the barycentric coordinates x Q to the 
black-hole coordinates x a . This transformation, between 
two systems of harmonic coordinates, can be fully worked 
out (Sec. V), relying on previous work by Kopcikin [19|, 
Brumberg and Kopcikin [20 , Damour, Soffcl, and Xu |9J, 
Kopeikin and Vlasov [21], and Racine and Flanagan [12] . 

The matching procedure determines the coordinate 
transformation completely (Sec. VI C), and it produces 
a justification of the earlier statement that the black hole 
can be treated as a post-Newtonian monopole. (A fuller 
discussion of this point is presented at the end of Sec. IV.) 
This statement, therefore, is a strict consequence of the 
Einstein field equations, rather than an artificial assump- 
tion. At first post-Newtonian order, the gravitational 
field of a black hole is that of a pure monopole, and it 
would be inconsistent to endow the black hole with an 
additional multipole structure. 

The matching procedure produces also an equation of 
motion for the moving black hole. It reads (Sec. VI E) 



d a U c: 



4(d a U b cxt - d b U: xt )v b 



+4d t U: xt +(v 2 -AU CKt )d a U e: 



-v a (iv b d b U cxt + 3d t U cxt ) 



+0(c- 4 ), 



(1.7) 



1 Here and below, the word "matter" describes a number of dif- 
ferent situations. The matter could be a continuous fluid, so as 
to model an accretion disk or a galactic core. It could also cor- 
respond to a collection of N — 1 bodies with weak self-gravity, 
making the black hole a member of an ./V-body system. Or else 
the domain ® could exclude a number N— 1 of small regions that 
would each contain a condensed body such as a neutron star or a 
black hole. In this last case, the post-Newtonian domain would 
contain no matter at all, but we will nevertheless refer to the 
N — 1 excluded regions as "matter." 



in which v = dz/dt is the black hole's velocity vector 
in the barycentric frame, and a = dv/dt is its accel- 
eration. The external potentials U cx t, U® xt , and ^cxt 
are defined as in Eq. (1.6), and they are evaluated at 
x = z(t) after differentiation. Equation (1.7! applies 



to a black hole moving in any post-Newtonian environ- 
ment. When this environment consists of (N — 1) ex- 
ternal bodies, so that the black hole is a member of 



an TV-body system, Eq. (1.7 1 reduces to the standard 
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(Einstein-Infeld-Hoffman) post-Newtonian equations of 
motion. These are listed, for example, in Exercise 39.15 
of Misner, Thorne, and Wheeler [25]. In effect, Eq. Oh 



states that the black hole moves on a geodesic of the met- 
ric of Eqs. ( 1.3 1-( 1.5 1, m which the (singular) potentials 
U, U a , and W are replaced by the (smooth) external po- 
tentials created by the distribution of matter outside the 
black-hole neighborhood SB. 

It is interesting to compare the differences between 
our derivation of Eq. ( |1.7| and the approach followed by 
Racine and Flanagan (RF) [12J. First, the work of RF is 
concerned with arbitrarily structured bodies (with weak 
or strong internal gravity), while our own work is con- 
cerned specifically with a nonrotating black hole, which is 
necessarily treated as a post-Newtonian monopole. Our 
work, therefore, is a specialization of theirs. Second, RF 
define the frame (i, x a ), which they call the body-adapted 
frame, by (essentially) setting the body's intrinsic mass 
dipole moment to zero; this is (essentially) the piece of 
goo that behaves as x a /f 3 . This coordinate choice does 
not, in general, constrain the tidal dipole moment; this 
is (essentially) the piece of goo that grows as x a . In our 
work, the coordinates (t, x a ) are defined so as to elimi- 
nate all mass dipole moments (both intrinsic and tidal) 
from the metric. This is made possible by the fact that 



we are dealing here with a specific type of body — a non- 
rotating black hole — instead of a general body whose 
nature is characterized only by an infinite set of mul- 
tipole moments. Indeed, the work of Zerilli [23 shows 
that in vacuum, an even-parity dipole perturbation of 
the Schwarzschild metric can always be removed by a 
gauge transformation; it is this gauge choice that de- 
fines our own version of the body-adapted frame, and 



the metric of Eqs. (1.1 1 and (1.2 1 reflects the complete 



absence of dipole terms. Third, in RF, the equations of 
motion are obtained by exploiting the integral form of the 
momentum-conservation identities that come as a conse- 
quence of the Landau-Lifshitz formulation of the Einstein 
field equations [24 . In our approach, the equations of 
motion are obtained directly by matching the black-hole 
and post-Newtonian metrics, and the computations are 
considerably simpler. This advantage is intimately tied 
to our complete control over the dipole terms; a deriva- 
tion of the equations of motion involving matching only 
would not be possible without the ability to set both the 
intrinsic and tidal mass dipole moments to zero. 

Finally, the matching procedure produces expressions 
for the tidal moments (Sees. VI D and E). In the barycen- 
tric frame they are given by 



-*<!. b 



-d ab U cxt + 1 [-d {ab) V cxt + Av c (d ab U^ - d c{a U^f) - 4d t{a U^ t - 2(v 2 - U cxt )d ab U cxt 
+3v c v( a d b ) c U C xt + 2v( a d b ) t U ex t + 3<9( a [4xt<9b)[/ C xt) + 0(c~ 4 ) 



(1.8) 



and 



B ab = 2e pq{a d p b) (U^ t - v q U ext ) + 0(c~ 2 ), (1.9) 



in which the external potentials are evaluated at x = 
z(t) after differentiation. The brackets around indices 
indicate symmetrization, while angular brackets indicate 
an STF operation: For any tensor A ab we have A/ ab -\ — 



\8„hA, where A 



\(A ab + A ba ) and A {ab) := A {ab ) 
S ab A ab . The expressions of Eqs. (1.8 1 and (1.9 1 are valid 



for any post-Newtonian environment. 

Given the vast difference in notations and ways of ex- 
pressing our results, we have not attempted to compare 
Eqs. ( |1.8[ ) and ( |1.9 | to the results obtained by Racine 
and Flanagan [12 , nor to those of Damour, Soffcl, and 
Xu [U Uni E]. We ca n sta te, h owever, that the spe- 
cialization of Eqs. (1.8 1 and (1.9 1 to an iV-body system 



is in perfect agreement with the corresponding results of 
Damour, Soffel, and Xu — see, in particular, Eqs. (4.29)- 
(4.31) of Ref. [TT]. We shall provide evidence for this 
statement in Sec. VII A. 

When the black hole is part of a binary system in cir- 
cular motion, the nonvanishing components of the tidal 



moments are given by (Sec. VII C) 
Gm! 



£ll-£: 



22 



b 3 

+ 0(c- 4 ) 
3GW 



1 + Tf 1 a(*Wc)' 

2(m + ml) 



(1.10) 



b 3 

+ 0(c- 4 ) 



3m + Am' . , 2 
1 - -. ^(Wc) 2 



£12 = 

Bia = 
B 23 = 



3Gr 



2b 3 
+ 0{c- A ) 
3Gm' 



2(m + m') 

cos2Lut, (1.11) 

3m + 4m' 2 
2(m + ml) 



b 3 
3Gn, 



sin2wi, (1.12) 
v rc iCosujt + 0(c" 2 ), (1.13) 



b 3 



-v le isinujt + 0(c 2 ), (1.14) 



where m is the mass of the black hole, ml the mass of 
the companion, and b := \z — z'\ the orbital separation 
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between the two bodies (in barycentric harmonic coor- 
dinates); the components £13, £23, Bn, B22, and B12 all 
vanish for circular orbits, and £33 = — (En + 822)- In 
these equations i> re i := \v — v'\ = \J G(m + to') jb stands 
for the (Newtonian) orbital velocity of the relative orbit, 
and 



G(m + to') 



b 3 



+ 0(c" 4 ) 



3m 2 + 5mm' + 3m' 2 2 

1 — : — (Wrel/C) 



2(m + m') 



A2 



(1.15) 



is the orbit's post-Newtonian angular velocity. Equations 
( 1.8 )-( 1.15) are expressed in the barycentric coordinates 



x a . In the main text (see Sec. VII D) we also list expres- 
sions that are valid in the black hole's moving frame. In 
the comoving coordinates x a we find that Eqs. (1.101- 



(1.141 are unchanged, except for the fact that the tidal 



moments must now be expressed in terms of the trans- 
formed phase variable ui. The transformed angular ve- 
locity Cu, given by 



G(m + to') 



b 3 

+ 0(c~ A ) 



3to 2 + 7mm' + 3to' 2 2 

-(Wrol/C) 



2(to + to') 5 



(1.16) 



accounts for the change in time coordinate as well as the 
geodetic precession of the moving frame relative to the 
barycentric frame. Notice the different coefficient in front 
of mm' in the post-Newtonian term. 

Equations (1.7l-(1.16l are the main results of this 
work. As an application we examine (Sec. VIII), in a 
suitable choice of gauge, the intrinsic geometry of the 
event horizon of a tidally-deformed black hole. We also 
calculate (Sec. IX A) the rate at which the black hole 
acquires mass as a result of its tidal interaction with the 
companion body. We find that this rate of tidal heating 
is given by 



Gto 



G(m+m') 



5c 15 (to + to') 8 
5to 2 + 12mm' 



b 

6m' 2 



1 



(to + to') 



r V 2 



(Urel/c) S 



+ 0(c- 4 ) 



(1.17) 



The rate at which the tidal coupling increases the black 
hole's angular momentum J can next be obtained from 
the rigid-rotation relation toc 2 = ujJ. 



C. Relevance of this work 

Our interest in this paper is mostly in the exploration 
of the tidal dynamics of black holes, in a weak-field con- 
text in which the tidal fields can be determined explicitly. 
These dynamics produce the tidal heating of the black 



hole, an increase in mass (and angular momentum and 
surface area) that is produced entirely by an influx of 
gravitational energy across the horizon. This fascinating 
phenomenon was studied before, most notably by Poisson 
and Sasaki [35], Alvi [SB], Price and Whclan [37], Hughes 
[28] . Martel [29], and Poisson [HI [30]. We provide here 
some additional insights. 

The tidal heating of a nonrotating black hole is gener- 
ally very small. Relative to the energy radiated away by 
gravitational waves, the effect is of order (v/c) 8 . In prac- 
tical terms, the effect is likely to be too small to be ob- 
served in a gravitational-wave signal that would be mea- 
sured by ground-based detectors such as LIGO, VIRGO, 
and GEO600. For example, Alvi [2B] calculated that for 
binary systems involving black holes with masses ranging 
from 5 to 50 solar masses, tidal heating is negligible: It 
contributes only a small fraction of a wave cycle during 
the signal's sweep through the detector's frequency band. 

In some circumstances, however, the tidal heating is 
a significant effect that should not be neglected [57] . In 
particular, it is likely to be observed in gravitational- 
wave signals that would be measured by a space-based 
detector such as LISA. For example, Martel [29] showed 
that during a close encounter between a massive black 
hole and a compact body of much smaller mass, up to 
approximately five percent of the total radiated energy is 
absorbed by the black hole, the rest making its way out 
to infinity. Hughes [28] calculated that when the massive 
black hole is rapidly rotating, tidal heating slows down 
the inspiral of the orbiting body, thereby increasing the 
duration of the gravitational-wave signal. For example, a 
1 Mq compact body on a slightly inclined, circular orbit 
around a 10 6 Mq black hole of near-maximum spin will 
spend approximately two years in the LISA frequency 
band before its final plunge into the black hole; Hughes 
shows that tidal heating contributes approximately 20 
days (and 10 4 wave cycles) to these two years. 

Another situation in which the tidal heating might be 
"measured" is in the numerical simulations of black- hole 
mergers. To reveal this small effect the simulations will 
need to be performed with very high accuracy and reso- 
lution, but this should become possible within the next 
few years. The simulations would reveal a steady growth 
in the irreducible mass of each black hole, at a rate that 



should be compatible with Eq. (1.171 if the holes are non- 



rotating. A recent paper by Boyle et al. [3T] suggests that 
the tidal heating of black holes might already have been 
seen in numerical simulations — witness their Figure 4. 

As an additional application of our work we mention 
another connection with numerical relativity. If the sim- 
ulations of black-hole mergers are to describe realistic 
situations of astrophysical interest, it is imperative that 
they proceed from data that correctly describe the ini- 
tial state of the system (the correct state of motion, 
and the correct amount of initial radiation). The idea 
that well-controlled initial-data sets could be constructed 
with the help of post-Newtonian theory is an old one 
[31 [331 [Ml 135] IM1 1ST] , but its initial formulation did not 
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account for the strong internal gravity of each black hole, 
which cannot be approximated by a post-Newtonian se- 
ries. Alvi [TnjES] was the first to remedy this situation by 
matching the post-Newtonian metric to black-hole met- 
rics, in the way that was reviewed in Sees. I A and I B. 
Alvi's work was improved upon by Yunes et al. |17j (see 
Sec. I A), and our work contributes an additional im- 
provement. While Alvi and Yunes calculated the black- 
hole tidal moments for circular motion only, and only 
to leading order in the post-Newtonian expansion, here 
we calculate £ a b to first post-Newtonian order, and B a b 
to leading order, for more general situations. Our re- 
sults could be used to generate an improved version of 
the Alvi- Yunes metric, which could then be used to con- 
struct improved initial-data sets for numerical relativity. 



D. Organization of this paper 

The technical portion of the paper begins in Sec. II 
with an introduction to the description of tidal environ- 
ments in terms of STF moments £ a b and B a b- In Sec. Ill 



we present the black-hole metric of Eqs. (1.1 1 and (1.2 1 



In Sec. IV we introduce the post-Newtonian metric of 



Eqs. (1.3 1— (1.5), and discuss the decomposition of the 



gravitational potentials int o a black-hole piece and an 
external piece, as in Eq. (1.6 1. In Sec. V we review 



the coordinate transformation between the barycentric 
frame (t, x a ) and the black-hole frame (t, x a ), and we cal- 
culate how the post-Newtonian potentials change under 
this transformation. The end result of this computation 
is a post-Newtonian metric expressed in the same coor- 
dinates as the black-hole metric of Sec. III. In Sec. VI we 
match the black-hole and post-Newtonian metrics and 
derive our expression of Eq. (1.7 1 for the black hole's ac- 



celeration vector, as well as Eqs. (1.8 1 and (1.9 1 for the 
tidal moments. In Sec. VII we specialize our results to 
the specific case of a two-body system. We first calcu- 
late the tidal moments for generic orbital motion, and we 
next specialize these results to circular motion; this gives 
rise to Eqs. (1.10)— (1.161 above. In Sec. VIII we present 



an application of our results: We examine the intrinsic 
geometry of the event horizon of a tidally-deformed black 
hole, in a suitable choice of gauge. And finally, in Sec. IX 
we apply our results to a calculation of the tidal heating 
of a black hole by an external body on a circular orbit; 
this is a gauge-independent effect. We first consider the 
case of a nonrotating black hole and obtain Eq. (1.171. 



We next consider the case of a rotating black hole; the 
result was not displayed above, but it can be found in 
Eq. (IO) below. 



Throughout the paper (except in Sees. Ill B and VIII) 
we work in quasi-Lorentzian coordinates x a — (ct, x a ) or 
x a = (ci, x a ), and we adopt a standard three-dimensional 
notation when we deal with spatial components. For ex- 
ample, we use v = (v x ,v v ,v z ) to denote a Cartesian 
vector with components v a in a flat, three-dimensional 
space. Indices on v a are manipulated with the Euclidean 



metric 5 a b, and e a bc is the familiar permutation symbol. 
The Euclidean norm of v is \v\ :— \J 5 a bV a v b , and we let 
r := | a; | and f := \x\. Because the paper is devoted to 
a post-Newtonian treatment of tidal gravity, we find it 
useful to use conventional units in which G and c are not 
set equal to one (Sec. Ill B is again an exception in this 
regard); as in much of the literature on post-Newtonian 
theory, we use c -2 as a formal expansion parameter. 



II. TIDAL SCALES AND TIDAL MOMENTS 

As shown in Sec. I B, the black hole's tidal environ- 
ment is described by the symmetric-tracefree (STF) ten- 
sors £ a b(t) and B &(i), and it is characterized by the 
length scale 1Z. Our purpose in this section is to for- 
mally introduce these quantities, and to set the stage to- 
ward the computation of the black-hole metric in Sec. III. 
The gravito-electric tidal moments £ a b can be introduced 
most simply in the context of Newtonian gravity; we shall 
do this first. The gravito-magnetic moments B a b do not 
exist in the Newtonian theory, but they appear in a rcla- 
tivistic description of tidal fields. (The overbar, we recall, 
indicates that the tidal moments are evaluated in a ref- 
erence frame that moves with the black hole.) 

In Newtonian gravity the total potential U can be ex- 
pressed as in Eq. ( 1.6 1, with the first term Gm/f describ- 



ing the black hole (here 



5 a bX a x ), and the second 



term ?7 ext describing the gravitational influence of the 
external matter. Assuming that the external potential 
varies slowly in the black-hole neighborhood SB, we ex- 
press it as a Taylor expansion in powers of f: 

U ext (i, x a ) = U cxt (t, 0) + g a (t)x a + l -£ ab {t)x a x b 

+ ^B abc (i)x a x b x c + 0(f A /TZ i ). 

Here g a := d a U cx t(i, 0) is the gravitational force (per unit 
mass) acting on the black hole, B a b := d ao U cx t(t,0) is 
the quadrupole moment of the external potential, and 
£abc '■= d a b S U ex t(i,0) is the octupole moment. The first 
term in the expansion does not depend on the spatial 
coordinates and plays no role in the gravitational inter- 
action of the black hole with the external distribution of 
matter. The second term g a x a also plays no role, be- 
cause we are working in a noninertial frame attached to 
the moving black hole. What remains is the pure tidal 
potential, t/ tidal = \£ ah x a x h + \£ abc x a x b x c + 0(f 4 /K 4 ), 
expressed in terms of the spatial coordinates and the tidal 
moments (which depend on time). Notice that the tidal 
moments are defined as fully symmetric tensors. And be- 
cause U ext satisfies Laplace's equation in SB (recall that 
the black-hole neighborhood is assumed to be empty of 
matter), the tidal moments are also tracefree. The quan- 
tities £ a b and £ a bc (and all higher-order moments) are 
therefore STF tensors: B a b = Bi ao ) and B a bc — £(abc)- I n 
addition, because c~ 2 U cxt is dimensionless, c~ 2 B a b has a 
dimension of inverse length squared, and the tidal length 
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scale TZ is defined such that the components of c~ 2 £ ab 
are typically of order TZ" 2 . In the preceding expression 
for Z7 ex t we have assumed that the components of c~ 2 £ abc 
are of order TZ~ 3 , and that each additional term in the 
expansion comes with an additional power of f /TZ. 

In general relativity the definition of the tidal moments 
£ab, £aba Bab, Babe (and higher-order moments) requires 
more refinement. The relevant tools were introduced by 
Thornc and Hartle [6] and Zhang [38] . We consider a 
neighborhood of a geodesic world line 7 in an arbitrary 
spacetime. In this neighborhood the metric is assumed 
to satisfy the vacuum field equations. (In Sec. Ill a black 
hole will be placed on this geodesic, and our spacetime 
will become the "external spacetime" of Sec. I A. For the 
time being, however, the black hole is absent.) We install 
a normal coordinate system (i,x a ) in the neighborhood. 
It possesses the following properties: (1) the spatial co- 
ordinates x a vanish on 7, and i is proper time on the 
geodesic; (2) the metric takes on Minkowski values on 
7; (3) all Christoffel symbols vanish on 7; and (4) the 
coordinates are harmonic, in the sense that the metric 
satisfies the conditions dp(^/—gg al3 ) — everywhere in 
the neighborhood. 

The metric near 7 admits an expansion in powers of f. 
By virtue of property (2) the zeroth-order terms are con- 
stant, and by virtue of property (3) there are no terms at 
first order. The terms at second and higher orders con- 
tain information about the curvature of spacetime near 
the geodesic, and it is those terms that describe the tidal 
environment around 7. This environment is character- 
ized by the scaling quantities TZ, C, and T, with TZ de- 
noting the radius of curvature on 7, £ the scale of spatial 
inhomogeneity, and T the time scale over which changes 
occur in the environment. The tidal environment is de- 
scribed precisely by the tidal moments £ a b(t), £abc(t), ■ • • 
and B ab (t),B abc (t),- ' ■, which are STF tensors that de- 
pend on t only. They are related to components of the 
Riemann tensor and its derivatives evaluated on 7 - 
see Zhang's Eqs. (1.3) for definitions [38 . The gravito- 
clcctric moments scale as 



(2.1) 



in which an overdot indicates differentiation with respect 



C 2 £ab r 


- nr 


2 


C &abc r 




2 L~ 


C~ 2 £ab r 




2 T - 



to t. It is the relations of Eq. (2.1 1 that define the scales 



TZ, C, and T. The gravito-magnetic moments scale as 



,-3 f 
-3» 



c 3 B ab ~ (v/c)TZ 



-2 

-2/.-1 



c s B abc ~ {v/c)TZ- z C- 
' 3 B ab ~ (v/c)TZ- 2 T-\ 



(2.2) 



in which v ~ C/T is a velocity scale. In a slow- motion 
context we have that v/c -C 1 , and the gravito-magnetic 
moments are smaller than the gravito-electric moments 
by a small factor of order v/c. 



To illustrate the meaning of these tidal scales we return 
to the example presented in Sec. I A, in which the tidal 
environment is provided by an external body of mass ml . 
We suppose that the geodesic is at a distance b from the 
external body. In this situation v/c ~ \J M' /b, where 
M' := Gm! /c 2 is the characteristic gravitational radius 
of the external body. The tidal scales are then given by 



TZ ~ yJb 3 /M', C ~ b, T ~ y/P/Gm'. 

We notice that C ~ (v/c)TZ and T ~ TZ/c; in a slow- 
motion situation we have that C -C TZ. 

The metric near 7 takes the form derived by Zhang 



.950 

90a 
9ab 



1 n£ab% £ a b c X XX -(-••-, (2.3) 

C z 6C Z 



2 
3c 



■p -b-c-d 



10 



&b c X X T & a foX 



(2.4) 



$ab ( 1 H - 9 £cd% ^"09 ^cde-E % 



3c 2 



-p 



1 



21c 4 \ ( a€b ) c P^ d% % — 5^ Cc P( a ^ 



(2.5) 



where e a bc is the permutation symbol. The neglected 
terms involve higher powers of f, and higher-order tidal 
moments. 



The first tidal term on the right-hand side of Eq. (2.3 1 
is of order (f/TZ) 2 , and the second term is smaller than 
this by a factor of order r/C; the neglected terms are 
smaller still, by additional factors of order r/C In 



Eq. (2.4 1 the first term is of order (v/c)(r /TZ) 2 , and the 
second term is smaller than this by a factor of order 
r/C; taking into account the scalings C ~ (v/c)TZ and 
T ~ TZ /c, t he same is true of the third and fourth terms. 
In Eq. (2.5 1 the first tidal term is of order (f/TZ) 2 , and the 



second term is smaller by a factor of order f/C; the third 
and fourth terms are smaller than this by another fac- 
tor of order f/C, and they come also with an additional 
factor of order (v/c) 2 . 

These considerations lead us to the following conclu- 
sion: If we restrict the neighborhood of 7 to be such 
that f is everywhere much smaller than C, then Zhang's 
metric can be simplified to 



,9oo 



9ab 



1 tj £ab% % 
C l 



O 



3c 3 



e abp B p c x b x c + O 



TZ 2 C)' 

v f 



3 \ 



: TZ 2 C J ' 



Sab[l+ \£cdX C X d )+0' ' 



TZ 2 C 



(2.6) 
(2.7) 
(2.8) 



This simplified form, which involves the lowest-order 
tidal moments only, and which neglects their time deriva- 
tives, shall be sufficient for our purposes below. 
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III. TIDALLY DEFORMED BLACK HOLE 



Metric of a deformed black hole 



The metric of Eqs. (2.6 1— ( 2.8 1 describes the tidal envi 



ronmcnt around a geodesic 7 in an arbitrary spacctimc, 



with the only restriction that the geodesic's neighbor- 
hood must be empty of matter. We now place a nonro- 
tating black hole of mass m on this geodesic, and modify 
the metric to account for its gravitational effects. As we 
shall show below (in Sec. Ill B), the metric of the tidally 
deformed black hole is given by 



500 

90a 
9ab 



r 2 (l -M/rfE* + 



K 2 £ 

2 -f\l - M/f)(l + M/f) 2 B a + o( V -^ 



1 - M/f 
1 + M/f 

1 4- M/f 

1 J,- n a n b + (i + M/f) 2 lab - r 2 (i + M/r) 2 s«n a n b 

1 — M I r 

- Mr (1 + M/f) 2 (l + M 2 /3f 2 )(n a £* + e*Sl b ) - f 2 (l - M/f) 2 (l + M/ff lab £* 

2^q 



(3.1) 

(3.2) 



Mf(l + M/f) 2 (l - M 2 /3f 2 )£2 b + O 



n 2 c 



(3.3) 



in the hole's neighborhood which is formally defined 
by f < r max with r max -C C (this is a refinement of the 
definition presented in Sec. I B). Here M := Gm/c 2 is 
the black hole's gravitational radius, il a := x a /f is a unit 
radial vector, and j a b '■= Sab — Qa^b- The tidal potentials 
are given by 



pq 

£ q 

°ab 



1 



£ c d^l c fl a 



cr 

^alb^ 



lab£ q 



1 



(3.4) 
(3.5) 
(3.6) 
(3.7) 



the label q stands for "quadrupole." The metric of 
Eqs. (3.1 1— ( 3.3 ) is presented in harmonic coordinates, 



and it is an approximate solution to the vacuum field 
equations linearized about the Schwar zsch il d me tric. In- 
deed, setting £ ab — B ab = in Eqs. (3.1 1— ( 3.3 > returns 
the Schwarzschild metric in harmonic coordinates, and 
the tidal potentials represent a pure-quadrupolc metric 
perturbation. It is easy to see that when f is much larger 
than M (but still much smaller than £), Eqs. ( 3. 1 1— ( 3.3 ) 



reduce to Eqs. ( 2 . 6 1— ( 2 . 8 1 . This shows that Zhang's met- 
ric provides the appropriate asymptotic conditions for the 
metric of a tidally deformed black hole; these replace the 
asymptotically-flat conditions that would be appropriate 
for an isolated black hole. Like the metric of Eqs. (2.6|- 
(2.8), the black- hole metric neglects tidal terms that are 



smaller than the dominant ones by additional factors of 
order f/C <C 1; within this approximation the metric is 
accurate to all orders in M/f. 



The method by which the metric of Eqs. (13. lb ( 3.3 > is 



obtained is explained in the next subsection. The reader 
not interested in those details can immediately proceed to 
Sec. IV, in which we present the post-Newtonian metric 
to which the black-hole metric will be matched. Before 
we proceed, however, it is useful to note that the com- 
ponents of the black-hole metric that are required for 
matching are 



1 ~ M/f 1 
9m - ~l + Af/f~^ (1 



M/f) 2 £ ab (t)x a x b 



O 



f 3 \ 



(3.8) 



and 



So a 



(1 - M/f)(l + M/f) 2 e abp B p c (t)x b x c 



3c 3 



■of-— Y 

1 n 2 L 



(3.9) 



We recall that M = Gm/c 2 is the black hole's gravita- 
tional radius. 



B. Derivation 

We begin this subsection with a warning on a change 
of notation; this change applies to this subsection only. 
To simplify the notation we shall refrain from displaying 
the over bar, in spite of the fact that we continue to work 
in a reference frame that moves with the black hole. In 
addition, we shall use relativist's units, in which c and G 
are both set equal to unity. And finally, in this section we 
use g a f3 to denote the Schwarzschild metric (as opposed 
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to the perturbed metric displayed in Sec. Ill A), and 
we let h a @ denote the tidal perturbation; the perturbed 



metric is therefore g a @ 

To obtain Eqs. pjV(|3J| we rely on Ref. [IS] 



which the metric of a tidally deformed, nonrotating black 
hole is presented in light-cone coordinates (v 7 p,8,(f>). 
The advanced-time coordinate v is constant on past 
light cones that converge toward p = 0, p is the usual 
Schwarzschild radial coordinate that measures the area 
of closed surfaces of constant (v,p), and 6 A = (6,(f>) are 
angular coordinates on these surfaces. In the light-cone 
coordinates the Schwarzschild metric takes the form 

g a/3 dx a dx f3 = -/ dv 2 + 2 dvdp+ p 2 fl AB dO A d6 B , (3.10) 

with / := 1 - 2M/p and Q AB d6 A dQ B := d0 2 + sin 2 6 dcf> 2 
denoting the metric on the unit two-sphere. The nonva- 
nishing components of the tidal perturbation are 



hi 



2 



,4 M = -yf(£ A ~B A) , 



(3.11) 
(3.12) 



/light 
AB 



1 



P \l-2M 2 /p 2 )£\ B 



+ -p 4 (l-6M 2 /p 2 )B AB , (3.13) 



in which the label "light" indicates that the perturbation 
is presented in the light-cone gauge. The tidal potentials 
are given in terms of the tidal moments £ a b{v) and B ab (v) 
by the relations 



£ q = s ab n a n\ 

£ A = fl A £ ab ^l b , 



AB 

B\ 



2Q a Sl R £ a \ } + fl AB £ q , 



^ A ^apq^>^ B ^1 



(3.14) 
(3.15) 
(3.16) 
(3.17) 



b\ b = n\e apq WB\n h B + n B e bpq WB\n% (3.18) 

where Sl a :— (sin cos cf>, sin 8 sin </>, cos 6) and Sl a A :— 
d£l a /39 A . The perturbed metric is g a p + h a p, and it is 
straightforward to show that this metric satisfies the vac- 
uum field equations linearized about the Schwarzschild 
solution. The metric is accurate up to terms of order 
(r 3 /1Z 2 £) that come from higher-order tidal moments, 
and terms of order (r 3 /1Z 2 T) that come from the time 
derivative of the quadrupole moments. For the purposes 
of substitution into the field equations, £ a b and B a b can 
be considered to be time-independent. 

The black-hole metric is expressed in the light-cone co- 
ordinates (v, p, 9 A ), and we wish to perform a transforma- 
tion to harmonic coordinates (t, x a ). We accomplish this 
in two steps. First, we perform a gauge transformation to 
change the perturbation h a p from its current light-cone 
gauge to a harmonic gauge. Second, we perform a back- 
ground coordinate transformation from the (background) 
light-cone coordinates to the (background) harmonic co- 
ordinates. This strategy, and its implementation detailed 
below, was suggested to us by Steve Detweiler (personal 
communication) . 



Let a set of four scalar fields be defined by 



T 
X 
Y 

Z 



v- p-2Mlog(p/2M - 1), (3.19) 

(p- A/) sin 6» cos 0, (3.20) 

(p- M) sin sine/), (3.21) 

(p-M)cos9, (3.22) 



and let us collectively denote the members of this set by 
X^\ It is straightforward to show that each one of the 
scalar fields satisfies the wave equation 



1 
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df,( g pPd a XM)=0 



(3.23) 

in the Schwarzschild spacetime. Here V Q is the covariant- 
derivative operator compatible with the Schwarzschild 
metric g a p, and g is the metric determinant. The state- 
ment of Eq. (3.231 is coordinate independent. When, 
however, we choos e (t — T, x = X y = Y, z = Z) as coor- 



dinates, then Eq. (3.231 becomes 9/3 (V — 99 ) — 0j the 



familiar statement of the harmonic coordinate condition. 



Equation (3.23) therefore provides a coordinate-invariant 
way of stating that the scalar fields X^' form a set of 
harmonic coordinates for the Schwarzschild spacetime. 

We now demand that X^' be harmonic coordinates 
for the perturbed spacetime, in addition to being har- 
monic coordinates for the Schwarzschild spacetime. To 



achieve this we rewrite Eq. (3.231 in terms of the full 
metric g a p + h a p and its associated covariant-derivative 
operator. We write g a P — h a P for the inverse metric and 
yj—g(\ + |/i) for the metric determinant, and we ma- 
nipulate the indices with the background metric; for ex- 
ample h a/3 = g ai g^ s h 1 s and h = g a ^h a p. After some 
straightforward manipulations we find that in addition 
to Eq. (3.231, the scalar fields must also satisfy the set of 



equations 



V a (f^V^W)=0, 



(3.24) 



where ip a /3 = h a p — \g a ph is the "trace-reversed" met- 
ric perturbation. This equation can be interpreted as a 
gauge condition on h a p: The perturbation will be in the 
harmonic gauge if it satisfies the conditions of Eq. (3.241. 



This equation is coordinate invariant, and the harmonic- 
gauge condition can be imposed even when the coordi- 
nates x a do not coincide with the scalar fields X^'. 



In Eq. (3.11 ) — (3.13 1 the metric perturbation is pre- 



sented in the light-cone gauge. A transformation to the 
harmonic gauge will be generated by the vector field £ Q , 
so that 



(3.25) 



It is straightforward to show that if is to satisfy 



Eq. (3.241, then the vector field must satisfy the four 



equations 



2(V Q ^ /3 )V Q V /3 A^ ) ; (3.26) 
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in our case the perturbation is traceless in the light-cone 



gauge, so V Q | 
vector 



1 



e P = oP 3 £ q , 



5 J-2 



p b f 



3(p-M) 



It may be verified that the 

(3.27) 
(3.28) 

£\+\p 2 {p 2 -&M 2 )B\ (3.29) 



is a solution to Eqs. 
tidal moments £ a b and B ab can be considered to be time- 
independent, because the time derivative of £ Q is smaller 
than its spatial derivatives by a small factor of order p/T. 
Sub stitution of Eq s. (|3.11|"(|3.13| and Eqs. \Z.27\ - 



(3.291 into Eq. (3.251 returns the tidal perturbation in 



the desired harmonic gauge. We obtain 



7^ harm 

i^ harm 
a vp 

i^ harm 
,l PP 



i harm 



k harr 



i harm 
11 AB 



-p 2 ! 2 £\ 
P 2 f£\ 

-2 P 2 £\ 
Mp 2 



3( P - My 



:(3p 2 



6Mp + ±M 2 )£\ 



3 r 
P 5 f 2 



(3.30) 
(3.31) 
(3.32) 

(3.33) 



(3.34) 



P 



M 
Mp 2 
3(/9 - M) 



(3p 2 - 6Mp+ 2M 2 )£\ B . (3.35) 



The full metric is next obtained by adding h a to g a p as 
given by Eq. (13.10 1. 



The tidal perturbation is now correctly expressed in 
the harmonic gauge, but the metric is still written in 
terms of the original coordinates (v, p,0,cj>). Our final 
step is therefore to perform a coordinate transformation 
from these coordinates to the harmonic coordinates (t = 
T,x = X,y = Y,z = Z). We carry this out in two 
stages. First, we effect a transformation from (y, p) to 
(t, r), leaving the angular coordinates alone; here t = 
v — p — 2A/log(p/2M — 1) is harmonic time and r = 
p — M is the harmonic radial coordinate. This coordinate 
transformation brings the Schwarzschild metric to the 
new form 



g a p dx a dx 13 



dt 2 



■M/r 



dr 2 



1 - M/r 
~l + M/r~" ' I- M/r' 
+ (r + M) 2 n AB d6 A dO B , (3.36) 



and the tidal perturbation becomes 



flight 



-r 2 {\ -M/r) 2 £ q , 
-r 2 (l + M/r) 2 £ q , 



(3.37) 
(3.38) 



/light 

/light 
n rA 
/light 
" AB 



= -r 3 (l~M/r)(l + M/r) 2 B\, 



(3.39) 



-Mr 2 (l + M/r) 2 (l + M 2 /3r 2 )£ q A , (3.40) 



-r 4 (l - Af/r) 2 (l + M/r) 3 fl AB 



Mr 3 (l + M/r) 2 (l - M 2 /3r 2 )£ A 



AB- 



(3.41) 



In the second stage we go from the spherical coordi- 
nates (r, 9 A ) to the associated Cartesian coordinates 
x a = rfl a , that is, x — r sin 9 cos <f>, y = r sin 9 sin cj), and 
z = rcos9. The transformation matrix is dr/dx a = fi a 



3.26b. In this computation the an . d dd 



l Q A , in which fl a :— S a b^ 



ant 



fl A := fl AB S a f)Cl B . The transformation of g a p + h a p un- 
der this change of coordinates is easy to carry out, and 
the end result is the metric of Eqs. (3.1 1-(3.3 1. 



IV. POST-NEWTONIAN METRIC 

In this section we take our attention away from the 
black-hole neighborhood 33, and we take the global view 
of the post-Newtonian domain @, which was introduced 
in Sec. I B. Let us recall that & is spatially limited by 
a sphere of radius r ncar centered on the post-Newtonian 
barycenter, and that this sphere marks the boundary of 
the near zone. The domain excludes a sphere of radius 
^min centered on the black hole, inside which the hole's 
gravity is too strong to be adequately described by post- 
Newtonian theory. We recall also that there exists an 
overlap 6 between the black-hole neighborhood 33 and 
the post-Newtonian domain '2. This is described by 
^min < r < r max , in which f max <C C is the boundary 
of 33. We assume that there is no matter in G . 

The metric in 9 is expressed in the post-Newtonian 
form 



3oo 



90a 



-l + -U+-(*-U 2 ) + 0(c- 6 ), (4.1) 



ib — 



;U a + 0(c- 5 ), 



1 + ^U)5 ab + 0(c-*), 



(4.2) 
(4.3) 



in which U is a Newtonian potential, U a a vector poten- 
tial, and \& a post-Newtonian potential. The metric is 
presented in harmonic coordinates x a — (ct,x a ). These 
coordinates are centered on the post-Newtonian barycen- 
ter, which defines the origin of an incrtial reference frame; 
we recall that they are distinct from the harmonic coor- 
dinates x a = (ci, x a ) used in 33. 

We assume that there is a distribution of matter within 
and we make only two assumptions regarding its na- 
ture. First, we assume that its gravity is sufficiently weak 
to be adequately described by the post-Newtonian met- 
ric of Eqs. ( 4.1 1— ( 4.3 1. Second, we assume that there 



is no matter in so that the black hole moves in an 
empty region of spacetime. Otherwise the distribution 
of matter within @ is completely arbitrary. It could be 
a continuous fluid, so as to model an accretion disk or 
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a galactic core; or it could correspond to a collection of 
N — 1 bodies with weak self-gravity, making the black 
hole a member of an TV-body system; or else the domain 
@ could exclude a number N — 1 of small regions that 
would each contain a condensed body such as a neutron 
star or a black hole. Our considerations allow for this 
degree of generality. 

We wish to solve the Einstein field equations in the 
vacuum domain & . To write them down it is convenient 
to express the post-Newtonian potential as 



^ = i> 



i d 2 x 

2 dt 2 '' 



(4.4) 



in terms of two new potentials ip and X. The degeneracy 
is broken by the field equations, which are given by (see, 
for example, Ref. [T2"] ) 



V 2 (7 

y2 [/ a 

X7 2 X 



o, 
o, 
o, 

2U, 



(4.5) 
(4.6) 
(4.7) 
(4.8) 



where V 2 is the usual Laplacian operator of three- 
dimensional flat space. The potential X is commonly 
referred to as the post-Newtonian superpotential. The 
vacuum field equations are augmented by the harmonic 
coordinate condition 



d t U + d a U a = 0, 



(4.9) 



which takes the form of a gauge condition on the gravi- 
tational potentials. 

Each one of the field equations is linear, and a solution 
representing a black hole and an external distribution of 
matter can be obtained by linear superposition. We treat 
the black hole as a post-Newtonian monopole of mass m 
moving on a trajectory described by the position vector 
x = z(t); we let v := dz/dt be the black hole's velocity 
vector, and a := dv/dt is the acceleration vector. In 
this treatment the solutions to the field equations are 
expressed as 



U(t,x) = 


Gm 


+ U ext (t,x), 


(4.10) 


\x-z\ 


U a (t,x) = 


Gmv a 


+ U? xt (t,x), 


(4.11) 


\x-z\ 


ijj{t,x) = 


Gm/i 


+ tp ex t(t,x), 


(4.12) 


\x - z\ 


X(t,x) = 


Gm\x — 


z\ + X cxt (t,x), 


(4.13) 



quantity n(t) that appears in -0 is a post-Newtonian cor- 
rection to the mass parameter m; this cannot be deter- 
mined directly from the vacuum field equations. 2 

The second terms on the right-hand sides of 
Eqs. ( 4. 10 1— ( 4.13 1 are the potentials created by the mat- 
ter distribution external to & . In & they sep ar atel y sat- 
isfy the vacuum field equations of Eqs. (4.5 1— (4.8|, and 



they separately satisfy the harmonic-gauge condition of 
Eq. (pi, 



d t U c 



8 T7 a 

u a u exi 



0. 



(4.14) 



We assume that the external potentials are smooth func- 
tions of the coordinates in a neighborhood of x = z(t). 
From Eq. (4.4) we find that the post-Newtonian po- 



tential ^ is given by 



Gm 
2^ 



Gm 



Gm 

< 

2s 



where s := x — z(t), s := \s\, and \& e 



S + *oxt, 

i 2 (415) 

'cxt = V'oxt + 2®t -^ext- 

The foregoing equations provide necessary and suffi- 
cient information regarding the post-Newtonian environ- 
ment of the black hole. A remarkable and important as- 
pect of our discussion is that the black hole is treated as a 
post-Newtonian monopole. This feature requires a justi- 
fication. At the level of Eqs. (4.10 1(4.13 I the monopole 



nature of the black hole is introduced as an assumption. 
As we proceed with the matching of the post-Newtonian 
metric with the black-hole metric of Sec. Ill, however, we 
shall see that our assumption is the only one that is con- 
sistent with the given structure of the black-hole metric. 
If we had instead given an arbitrary multipole structure 
to our black hole, the matching procedure would eventu- 
ally force us to set all higher multipole moments to zero; 
the only surviving moment is the black-hole mass m. As 
we shall see, therefore, the statement that the black hole 
is a post- Newtonian monopole is a strict consequence of 
the field equations of general relativity; no other multi- 
pole structure is possible. The monopole structure is not 
assumption introduced for its simplicity; it is a direct 
outcome of the matching procedure. 



V. 



TRANSFORMATION FROM BARYCENTER 
FRAME TO BLACK-HOLE FRAME 



The post-Newtonian metric of Eqs. (4.1 
the potentials of Eqs. (|4.10[), (|4.11[), and I 



4.15) 



1431), with 
is not 



where \x — z(t)\ is the Euclidean distance between the 
field point x and the black hole. 

The first term on the right-hand side of each equation 
represents the black hole, and we see that indeed, each 
black-hole potential has a monopole structure. The fact 
that the vector v appears in the vector potential U a is 
a consequence of the gauge condition of Eq. (4.9 1. The 



2 If the black hole were treated as a point particle, solving the field 
equations in the presence of a distributional energy-momentum 
tensor would reveal that /j = |ti 2 — f/ C xt(t, x = z), where 
v 2 = \v\ 2 . This is indeed the correct expression, but in our 
approach the determination of /i will come at a later stage, from 
a careful matching of the post-Newtonian metric with the black- 
hole metric of Sec. III. For the time being fi(t) will remain un- 
determined. 
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yet ready to be matched to the black-hole metric of 
Eqs. (3.1 1 — ( 3.3 k While both metrics are defined in the 



overlap region 6 and describe the same physical sit- 
uation, they are expressed in different coordinate sys- 
tems: The post-Newtonian coordinates (t,x a ) are de- 
fined everywhere in 3$, and they are attached to the 
post-Newtonian barycenter; the black-hole coordinates 
(i, x a ) are defined only in 3), and they are attached to 
the moving black hole. Each coordinate system is har- 
monic, and since they overlap in G there exists a coordi- 
nate transformation between them. This transformation, 
from harmonic coordinates to harmonic coordinates, is 
worked out in Sec. V A, following closely the treatment 
of Racine and Flanagan [12] . Their work was built on the 
post-Newtonian theory of reference frames developed by 
Kopeikin, |19j . Brumberg and Kopeikin |20j . Damour, 
Soffel, and Xu [9], and Kopeikin and Vlasov [2"Tj. In 
Sec. V B we calculate how the potentials of Eqs. (4.101, 



(4.111, and (4.151 change under the transformation from 



the barycenter frame to the black-hole frame. And fi- 
nally, in Sec. V C we express the transformed potentials 
in terms of irreducible quantities that will facilitate the 
matching with the black-hole metric, to be carried out in 
Sec. VI. 



A. Coordinate transformation 

The most general coordinate transformation that pre- 
serves the post-Newtonian expansion of the metric is 
given by [T^] 



t = t 



\a(t,x a ) + \l3(i,x a ) + 0(c- 6 ), (5.1) 



1 



x a + z a {t) + ^h a {t,x a ) + 0(c- 4 ), 



where 

a{t,x a ) = A(t) + v a x a , 



(5.2) 



(5.3) 



h a (i,x a ) = H a {t) + H a b {t)x b + -H a bc (t)x b x c , (5.4) 



with 



H ab (t) 
HabciJ) 



£abcR c (t) + ^v a v b -5 ab (A-^v 2 y (5.5) 



(5.6) 



The functions A, z a , H a , and R a arc freely specifiable 
functions of time t, while /3 is a free function of all the 
coordinates; these functions characterize the coordinate 
transformation. An overdot indicates differentiation with 
respect to t, and we have introduced 

v a := z a , a a := v a = T . (5.7) 

As before indices are raised and lowered with the Eu- 
clidean metric S a b, and we let v 2 = S a bV a v b - 



The transformation of Eqs. (5.1 1 and fl5.2| preserves 



the post-Newtonian expansion of the metric, but it does 



not necessarily keep the coordinates harmonic. To pre- 
serve this also we must set 

(3(i,x a ) = l -A5 ab x a x h 



+ (Saba-c + 5 ac a b + S bc a a )x a x b x c 
+ l{i,x a ), 



(5.8) 

and 7 is required to satisfy Laplace's equation: V 2 7 = 0, 
with V 2 denoting the Laplacian operator in the coordi- 
nates x a . 

Under the coordinate transformation the potentials be- 
come 



U(t,x a ) = U - A+^v 2 - a a x a , 



where 



(5.9) 



U a (t, x a ) = U a - Uv a + - ( V a + V a b x b + -V a bc x b x c 



(5.10) 



^(i,x a ) = 4>-4:U a v a + 2v 2 U+(A + v b x b )dtU 

1 

2' 



(f» + F a b x b + '-F a bc x b x c )d- a U 
G + G a x a + —G a bX a x b + —G a b c x a x b x c 



dry, 



(5.11) 



V a = (2A - v 2 )v a - H a + e a bc v b R c , (5.12) 



V\ = ^v a a b +^a a v b + 6 a b (^A-2v c a ( 

a f>c 

e bc n , 

V\ c = \{6\a c + 8 a c a b )-~8 bc a a , 

pa = H a_ Av a^ 

F a h = -5 a b (A- l -v 2 )- l -v a v b + e\ c R\ (5.16) 

F a bc = -{5 a b a c + 6 a c a b )+a a 6 bc , (5.17) 

(5.18) 



(5.13) 
(5.14) 
(5.15) 



G = l -A 2 -Av 2 + l -v A + H a v a , 

G a = [A - ^ 2 ja Q . - (A - ^v c a^jv a 

-e abc v b R c , (5.19) 
Gab = a a a b - v a a b - a a v b + 5 ab (v c a°) 



U ab A^\ 



G a bc — — g(<5aba C + 8 ac CL b + Sbc'da) 



(5.20) 
(5.21) 



Here A (3) stands for d 3 A/dP. 

The "hatted" potentials are equal to the original po- 
tentials evaluated at time t — t and position x a = 
x a + z a (t). For example, 

U(t,x a ) :=U(t = t,x a = x a + z a (t)). (5.22) 
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Because U now possesses, in addition to its original ex- 
plicit time dependence, an implicit time dependence con- 
tained in z a (t), some care must be exercised when taking 
time derivatives of U. We have, for example, 



in which the vector s := x — z(t) has been identified 
with x; we continue to use the notation f :— \x\ = 

(<Wi h ) 1/2 - 

The external potentials are given, for example, by 



dU 



dU 

~dt 



dx a 



t—t,x—x-\-z 



(5.23) 



dU f dU 
dx a ~ \dx% t=t , x=s+z 

The harmonic coordinate condition reads 

d{U -v a d a U + d- a U a = (5.24) 

when it is expressed in terms of the hatted potentials. 



B. Post-Newtonian potentials in the black- hole 
frame 



Our task in this subsection is to transform the poten- 



tials of Eqs. (4.101, (4.11), and (4.151 from the barycenter 
coordinates (t, x a ) to coordinates (t, x a ) that are centered 
on the black hole. Each coordinate system is harmonic, 
and the transformation was described in the preceding 
subsection. The most important pieces of the coordinate 
transformation are the functions z a (t), and for these we 
choose 



z a (t) = z a (t = t), 



(5.25) 



where z(t) is the black hole's position vector in the 
barycentric frame. In words, the coordinate displace- 
ments z a (t) are given by the black-hole position in the 
barycentric system, evaluated at the time t = i. With 
this choice it follows that the quantities v a and a a that 
appear in the coordinate transformation are the same as 
those contained in the post-Newtonian potentials. The 
other freely specifiable pieces of the coordinate transfor- 
mation will be determined in due course. 

The transformed potentials U, U a , and ^ are expressed 
partly in terms of the hatted potentials. We have, for ex - 
ample, U(i, x a ) = U(t,x a +z a ). According to Eqs. |4.10|, 



(4.111, and (4.15), the hatted potentials are 



U 



Gm 
f 

Gmv a 
f 

Gm 



4- TI a 



(5.26) 
(5.27) 



-„-b Gn 

—rv a v b x x + — 
Zr a r 

Gm 
' 2F 



( 1 2\ 



■a a x 



U cxt (t, x a ) := U cxt (t = t, x a = x a + z a (t)) . (5.29) 
They satisfy the harmonic coordinate condition 

dtUext ~ V a d a U cxt + daU^ t = 0. (5.30) 

Because it is well-behaved near x a , each external poten- 
tial can be expressed as a Taylor expansion about x = 0. 
For example, 

U ext (i, x) = U cxt (t, 0)+d- a U cxt (t, 0) + ^d- ai U cxt (t, 0) + - ■ • . 

(5.31) 

This defines the strategy behind our calculation of the 
transformed potentials: Each quantity that is smooth at 
x = will be expressed as a Taylor expansion. The 
potentials, therefore, will contain a piece that is singular 
at x = 0, and a smooth piece that will be expressed as a 
Taylor series. All Taylor expansions will be truncated at 
quadratic order, and they will involve derivatives of the 
external potentials evaluated at x = 0. 

The harmonic function jit, x a ), introduced in 
Eq. (5.8 1, is smooth at x = 0, and it also can be ex- 



pressed as a Taylor expansion. We write 



7 (*",x a ) = C{t)+ la {t)x a + - lab {t)x a x h 



G 



rT\ —a—b—c 
Jabc{t)X X X ■ 



(5.32) 



To ensure that this is a solution to Laplace's equation, we 
must choose the expansion coefficients to be symmetric- 
tracefree (STF) tensors. We express this property as 
lab = 7<afc) an d labc = J{ a bc) ■ The expansion coefficients 
are otherwise arbitrary, and they will be determined in 
due course. 

After a lengthy computation we obtain 



U 
U a 
f 



Gm 



+ U+ 1 U a x a + 2 U ab x a x b + 



U a + iU a b x b + 2 U a bc x' J x- + 



(5.33) 
(5.34) 



Gm 

hn*H 



Av a )x a + + A-2v 2 ) 



r* a X a 



J^ab X a X b 



(5.35) 



(5.28) with 
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JJ 

JJ ab 
U a 

■2 U be 
2*ab 



— a a + ds,U ex t, 

i7ext - * a &*t + \{2A - v 2 ) - V + \e\ c v h R c + \r, 
d b U« xt ~ v a d- b U cxt + h a a b + ±a a v b + \t\(^A - 2v c a c ) - \ 
\d b - c U^ ~ \v a dbcU C xt + ^(6%a c + S a c a b ) - ±a a 5 bc + ^ 

*cxt - ^v a U^t + 2v 2 U cxt + AdJJ ext + (H a - Av a )d- a U cxt + X -A 2 - Av 2 + ^v 4 + 

/5 .\ * i 
dafcxt - 4v b d a U b xt + y-v 2 - A}d- a U cxt - -v a v b d- b U cxt + v a dJJ cxt + Ad m U cxt 

+ (H b - Av b )d- ab U cxt +(A- \v 2 )a a -[A- ^v c a c y a - e abc (d b U cxt R c + v b R c ) - %, (5.43) 



e a hr R c 



77b, 



(5.36) 
(5.37) 
(5.38) 

(5.39) 

(5.40) 

(5.41) 

H a v a - C, (5.42) 



ab^cxt 



2v cdabU ext 



+ (|« 2 - A^d ab U ext - l,v c v (a d- b)£ U cxt + v (a d- b) JJ cxt + ^Ad m U cxt 
- a( a <%)0ext + ^S ab a c dJJ cxt - e c p{a R p d b) JJ cxt + ^(H c - Av c )d sab U cxt + ^a a a b - f( a a & ) 
+ ^8 ab (v c d c ) - ^6 ab A (3) - l -^ ab . 



(5.44) 



It is understood that in these expressions, the external potentials are evaluated at x — after differentiation. This 
notational convention will be retained below. ^ 



C. Decomposition into irreducible pieces 

To facilitate the matching procedure it is useful to de- 
compose the tensors JJ ab , JJ ab , JJ abc , and j & ab into 
their irreducible components. Equation (5.381 reveals 



that 2 U ab is already a pure STF tensor, because V 2 t/ cxt 
( ) everywhere near the black hole. 
We write 

iU ab = - S ab JJ + JJ {ab) + JJ [ab] , (5.45) 

which is a decomposition of JJ ab into a trace part, an 
STF part, and an antisymmetric part. Equation (5.40) 
implies 



metric part, and we write U abc = U( abc ) + V abc , where 
Vabc = l(U[ab] c + U[ ac ] b ) is what is left over of U abc after 
complete symmetrization. The first term is decomposed 
into trace and STF parts. For the second term we note 
that ^U[ ab ] c possesses 3x3 = 9 independent compo- 
nents, so that it can be expressed as e abp X p c , in terms 
of a general 3x3 matrix X ab . This, in turn, can be 
decomposed as X ab = ^S ab X + V ab + X[ ab ], in terms of 
trace, STF, and antisymmetric components. Finally, we 
write X[ ab ] = e a bpV p , which relates the 3 independent 
components of X[ ab ] to those of a vector V a . Altogether, 
we have §E/[ a f,] c = \e a b c X + S ac V b - 6 bc V a + e abp V p c . This 
produces V abc = S ab V c + S ac V b ~2S bc V a + e abp V p +e acp V p b , 
and we obtain the decomposition 



JJ 

lU(ab) 
lU[ab] 



A - dJJ cxt - v c a c , (5.46) 
5 <a^6) xt - V( a d b) U cxt + 2 V (a a b) 

1 

-<9[ a C4 C ] Xt - V[ad- b] U cxt + ^v [a a b] 

1 



(5.47) 



;£abcR C 



(5.48) 



To decompose JJ abc we first isolate its completely sym- 



2U abc = 2 U ( abc ) + -(8 ab 2 U c + 8 ac2 U b + 8 bc2 U a) 

+ 8 a b 2 V C + S ac 2 V b - 2S bc 2 V a 
+ tabp 2 V p r + e a 



V p 

^acp a' b : 



where JJ a := S bc 2 U( abc y The 18 independent compo- 
nents of JJ abc have been packaged into the 7 components 
of 2 C^(a6c) > the 3 components of JJ a , the 3 components of 
2 V a , and the 5 components of 2 V ab . Calculation shows 
that 2 V a — \ 2 U a and we obtain, finally, 



1G 



with 



9 3 

lU abc = 2U ( abc ) + ^(<5a6 iV c + $ac 2U (,) — ^S bc 2^a + e ahp 2^ P C + 6 acp 2 ^ P fc , 
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(5.49) 



2 U a 
,V ab 



2U (abc) 



-6 £ (a P9 %(^ Xt --^e xt ), 
(a a6 C> c OXt + 9ac^6 Xt + dbfiT) ~ I (vcd-abUcxt + V b d aS U eKt + V a d bE U cxt ) 



6 

+ ^di(5abdcU C xt + S ac d b U cxt + dbcd a U cxt ) + gTabc- 



(5.50) 
(5.51) 

(5.52) 



That the right-hand side of Eq. ( 5.52 ) is STF follows from of Eq. (5.30 1; and (iii) 7 a & c is itself an STF tensor, 
the facts that: (i) the potentials U ex t and U^ x t satisfy 

Laplace's equation; (ii) they obey the harmonic condition The decomposition of $ ab is 



2^ab = 2*(ab) 



1 



(5.53) 



with 



2^(ab) 



l(dtiU C xt + a 2 + v a a a -A ( V), 

(Z 2 

\2 V 



2 



2v c d- a - b U^ t 



^)d ab U cxt - ^v c v {a d- b)s U cxt + v {a d b)i U cxt + )^Ad lab U c - 
a(adh)U cxt + ^<5 Qb a c 9 5 C/o X t - e c p{a R p d b) Jj cxt + \{H C - Av c )d sab U cxt - ^S ab duU cx t 



1 



1 



+ 2 a (a a b) - v (a a b) - ^lab- 



(5.54) 



(5.55) 



The right-hand side of Eq. ( 5.55 ) is STF by virtue of the 



field equation satisfied by ^ cx t . In the barycentric frame 
we have V 2 ^^ = 2dttU cy ±] transforming to the black 
hole's moving frame gives instead 



V 2 * 



• xt — 9 tt U cxt — a a d a U ex t 



2v a d T „U . 



+ v a v"d ab U ext . 



(5.56) 



VI. MATCHING THE BLACK-HOLE AND 
POST-NEWTONIAN METRICS 

Because the black-hole and post-Newtonian metrics 
are now expressed in the same coordinate system (t, x a ) 
in the overlap region ff, we are finally ready to compare 
their expressions. We recall that the black-hole metric 
is valid when r <C C, while the global post-Newtonian 
metric is valid when f 3> M; the expansions provided 



f <^ C The metrics can be compared directly when f is 
such that M <C r -C C. A typical value of f in G might be 
r c ~ \J ML, and we have that f c /L ~ M/f c ~ v/c\ each 
quantity is indeed small. We write the metrics in Sec. VI 
A in a form that is ready for matching, and the matching 
conditions are extracted in Sec. VI B. In Sec. VI C we 
use them to determine the free functions associated with 
the coordinate transformation; this procedure reveals, in 
particular, the equations of motion for the black hole. In 
Sec. VI D we determine the free functions that appear 
in the black- hole and post-Newtonian metrics; it is here 
that the tidal moments £ ab and B ab are finally obtained. 
In our last subsection, Sec. VI E, we transform the equa- 
tions of motion and tidal moments from the black-hole 
frame back to the barycenter frame, in which they are 
most easily interpreted. 



in Eqs. ( 5.33 H( 5.35) are also restricted to the domain 
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A. Metrics 

In the overlap region G the black-hole metric can be 
expressed as a post-Newtonian expansion. Writing 



^^ + 0(c-) + 0(4) 



(6.2) 



and 



1 



^ = 4 N b + ^4™ + o(c- 4 ), 



(6.1) 



Soa 



3c 3 



e abp B p c (i)x b x c + 0(c- 5 ) + 



cK 2 £ 



(6.3) 



we get 



.9oo = -1 + 



2Gm 2G 2 m 2 1 



^pN-a-f, _ J_oPN-a-b 



The post-Newtonian metric is obtained by inserting 
Eqs. < pr^ -( pT35} within Eqs. j4l|)-([43]| ■ The result 
is 



2Gm 2Gm - r . s „ 2Gm , ; „ 
.900 = -l + -5^ 1^ ff«-^ a x a + ^^(m + A-2w 2 ) 



2G 2 ™ 2 



and 



4Gto 
2 



( C/ + JJ a x a + 2 U ab x a x b ) + j> (V + i tf - ^( U) 2 ^j 



1 2 \ 2 

c* \ &■ & I c 

.3 



2 I aUa, 



b ' 5 



0(0 + O 



90a 



K 2 £ 



{oU a + JJ ab x h + 2 U abc x b x c ) + 0(c' 5 ) + O 



cK 2 £ 



(6.4) 



(6.5) 



Comparing Eq. (6.4 1 to Eq. (6.2), and Eq. (6.5) to and 



Eq. (6.3 1 produces a complete set of matching conditions. 



B. Matching conditions 



From the absence of ^"-independent terms in Eq. (6.2 ) 
we get U + c~ 2 [Jii — ( U) 2 ] = 0(c~ 4 ), which implies that 
U = 0(c~ 2 ). The equation simplifies to 



U- 



1 



* = O(c- 4 ). 



(6.6) 



From the absence of terms linear in x a we get JJ a 



l [$> a - 2 UJJ a ] = 0(c- 4 ), which implies 



1 



(6.7) 



Comparing the singular terms in <?qq and taking into ac- 
count the facts that U — 0(c~ 2 ) and JJ = 0(c~ 2 ), we 
obtain 



H a - Av a = 0(cT 2 ), 



H + A-2v 2 = 0(c~ 2 ), 



(6.8) 



(6.9) 



£ a b — ab \ 



(6.10) 



the last equation can be valid if and only if 2 V ab is an 
STF tensor, a property that was already established in 
Sec. V C. Equation (6.101 is recovered by matching the 



terms that are quadratic in x a in 550 > an d this also reveals 
that 



^PN _ r, fly 

tab — *2*a&- 



(6.11) 



We observe that unlike all preceding equations, 
Eqs. (6.10 1 and (6.11) do not include an error term 
0(c~ 2 ); this is because each member of these equations is 
defined so as to possess a specific post-Newtonian order. 



From the absence of ^"-independent terms in Eq. (6.3 1 
we get 

U a = 0(c~ 2 ), (6.12) 
and from the absence of linear terms we obtain 

1 U ab = 0(c- 2 ). (6.13) 
Finally, matching the quadratic terms in produces 
e pa ( b B p c) = -6 2 U abc + Q{c- 2 ). (6.14) 
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The matching conditions of Eqs. (6.6l-(6.14l allow the 
determination of the quantities 

A, Z a , H a , R a , C, 7a , 7o0 , 7 afcc 

that appear in the transformation from the barycentric 
frame (t, x a ) to the black hole's moving frame (i, x a ). 
They allow also the determination of the functions 

Mi c abi c ab ' °ab 

that appear in the post-Newtonian and black-hole met- 
rics. 



C. Determination of the coordinate transformation 



From Eq. (6.6) we learn that JJ = Q(c ), and with 



the expression for JJ given in Eq. (5.36|), we find that 
A(i) is determined by the differentia 



equation 



A= -v 2 

2 



(6.15) 



where U cxt = U ext (t,0) is the external Newtonian poten- 
tial evaluated at x a = 0. With A determined, Eq. (6 
implies 



H a = Av a + 0(c- 2 ), 



(6.16) 



and this dete rmin es H a (t). 

From Eq. (6.7l we learn that JJ a = Q( c~ 2 ), and with 
the expression for JJ a given in Eq. (5.371, we obtain an 



expression for the black hole's acceleration vector: 

a a = d a U cxt + 0(c- 2 ), (6.17) 

where the external potential is evaluated at x a — after 
differentiation. This is a Newtonian approximation to the 
acceleration vector, and the post-Newtonian corrections 
will be determined below. 

With JJ = 0(c~ 2 ) Eq. flaB implies * = Q(c~ 2 ), and 



taking into account Eqs. (6.15 1— ( 6 . IT I , Eq. ( |5.42| reveals 
that 



C 



+ A(dJJ cxt + v a d a U cxt ) + -v A + 0(c~ 2 ), (6.18) 



in which each external potential is evaluated at x a = 
after differentiation. This equation determines C(t). 

Equation (6.131 implies that each irreducible pi ece of 
JJ ab must vanish to order c . According to Eq. 
we must have JJ = = JJ ( a b) — JJ l ab]- W ith Eq 



5.45) 



5.46) 



we reproduce Eq. (6.151. From Eqs. (5.47) and ( 6.17| we 

get 



lab = -AdptJfi* + 2v {a d- b) U cxt . 



(6.19) 



And from Eq. ( 5.48 1 we find 

Cabc-R = — 49[ a C/£] Xt - 



3^[adb]^ext, 



an equation that determines R a (i). 

Tak ing into account Eqs. ( |6.15| -( [6Tl7| , Eqs. ( |6.12| and 
KM imply 



7a 



_ 4 ^ext + (A y 2 + 3 u ext y a + Ad a U cxt 
-e abc v b R c + 0(c- 2 ). (6.21) 
We may now determine the post-Newtonian correc- 



tions to the acceleration vector. We return to Eq. (6.7 1 



in which we insert Eqs. (5.371 and (5.431. We next in- 



corporate Eqs. ( 6. 15 1— ( 6.17 1 , as well as Eqs. (6.201 and 



(|6.2l|. After simplification we obtain 
5 a * ext - 4v b d a U b ex 



a a = gajje 



1 



+ (v 2 ~4U cxt )d a U cxt 
+ 0(c- 4 ), 



(v c d E U cxt + 3dtU cx t)v a 
(6.22) 



where (as always) the external potentials are evaluated at 
x a = after differentiation. Equation (6.221 is a system 



of second-order differential equations for the functions 
z a (t); they represent equations of motion for the black 
hole. 

The last piece of the coordinate transformation that 
must be determined is 7 aoc . The information comes from 
Eq. (6.141 and the decomposition of Eq. (5.49). Com- 



paring the equations reveals that JJ a and JJ ( abc ) mus t 
both vanish. The first statement reproduces Eq. (6.171, 
while the second implies 



labc 



+ I {vcd-aiU^ + v b d a5 U cxt + v a d- bs U ext ) 
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di(S ab d s U cxt + 8 ac d b U cxt + S bc d a U eKt )- 



This can be expressed more compactly as 

labc = -±{d {ab Uf - V { ad b5) U^). 



(6.23) 



The coordinate transformation is now completely deter- 
mined by the matching conditions. 



D. Determination of the metric functions 



Equations (6.9 1 and (6.151 imply 



Uoxt, 



(6.24) 



in which the external potential is evaluated at x a = 0. 
As was first pointed out in the footnote that follows 
Eq. (4.131, this assignment can also be obtained by calcu- 



lating the post-Newtonian potential ip sec Eq. (4.4 1 



(6.20) for a point particle of the same mass as the black hole 
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nate the trace part of S> ab , we obtain 



Our matching procedure shows that // keeps the same 
value when the particle is replaced by the black hole. 

The electric components of the tidal fields are de- 
termined by Eqs. (6.10) and (6.11 1. After inserting 

Eqs. (5.381 and (5.53) and invoking Eq. (6.151 to elimi- and £™ = —2j$!/ ab \. An explicit evaluation of this yields 



(6.25) 



cPN _ 
°ab — 



-d(ab)^ cxt 



4v c d ab U^ t - Ach t{ <fitf - 2(v 2 - U cxt )d ab U CKt + v c v {a d b)E U cxt + 2v {a ch m U c 

I 



(6.26) 



where, as usual, the external potentials are evaluated at 
x a = after differentiation. The complete tidal poten- 



tials are £ a b 



in Eq. (6.1) 



°ab 



c S^b + 0(c ), as was expressed 



The magnetic components of the tidal fields are deter- 
mined by Eq. (6.14) and the decomposition of Eq. (5.491 



Taking into account the facts that 2 U a and 2 U i abc ) must 
both vanish (as was noted previously), comparing the 



equa tions reveals that B a b 
Eq. this is 



-l2 2 V ab + 0{c' 2 ). With 



B ab = 2e^ q d- b)p (U^ - v q U CKt ) + 0(c~ 2 ). (6.27) 

The metric functions are now fully determined by the 
matching conditions. 



E. Transformation to the barycentric frame 

In the moving frame (t, x a ) the black hole is situated 
at x a = 0. According to Eqs. ( 5.1 1 — ( 5.4 1, the position of 



the black hole in the barycentric frame is described by 
the parametric equations 



1 



ibh = t+^A(t)+0( C - 4 ), 



(6.28) 



bh 



= z a (t) + ^H a (t) + 0{c- 4 ). 



The equations of motion for the black hole, expressed 



in the barycentric frame, are obtained from Eq. (6.22 1 by 
replacing the hatted potentials (U cxt and so on) with the 
original potentials (U cx t and so on) using the correspon- 
dence of Eqs. ( |5.22[ ) and fl5.23| . Noting that Eq. ( |6.22 l 
is to be evaluated at i = tbh = t, we get 



= d a U c 



1 



d a ^ t - 4(cW c b xt - d b U? xt )v b 



+ Ad t U^ t + (v 2 - 4U cxt )d a U cxt 



v a (Av b d b U e 



3d t U c: 



+0(c" 4 ). (6.30) 



The external potentials were introduced in Eqs. (4.101, 
(4.111, and (4.151, and here they are evaluated at x = 
z{t) after differentiation. Equation (6.301 applies to a 
black hole moving in any post-Newtonian environment. 
When this environment consists of (N — 1) external bod- 
ies, so that the black hole is a member of an iV-body sys- 
tem, Eq. (6.30 1 reduces to the standard (Einstein-Infcld- 
Hoffman) post-Newtonian equations of motion. Because 
this connection is well understood, we shall not provide 
here a derivation of this well-known fact; the EIH equa- 
tions are listed, for example, in Exercise 39.15 of Misner, 
Thorne, and Wheeler [22] . 

Following Racine and Flanagan |12j we define barycen- 
tric tidal moments £ a b(t) and B a b(t) that are related to 
those of the black-hole frame by the transformation 



The first equation can be approximately inverted as 
t = ibh — c~ 2 A(tbh) + 0(c -4 ), and substitution into the 



second equation yields x^ h = z a (t ^) + c [H a (thh) — 
A(t hh )v a (t hh )} + 0(c- 4 ). With Eq. (|6.16| this becomes 



zg h = z a (i bh ) + 0(c- 4 ), 



(6.29) 



which is the same statement as Eq. (5.251. The position 



of the black hole in the barycenter frame is therefore ob- 
tained simply by evaluating the functions z a (t) at the 
time i — tbh- From this observation it follows that the 
black hole's barycentric velocity is v a (tbh), and its ac- 
celeration is a°(ibh)- Henceforth we shall omit the label 
"bh" on the barycentric time coordinate. 



£ ab (t) := M a c (t)M b d (t)£ cd (t), 



B ab (t) := M a c (t)M b d (t)B cd (t), 



where 



M ab (t) := S ab + ^e abc R c (i) + 0(c~ 4 ) 



(6.31) 



(6.32) 



is a post-Newtonian rotation matrix that accounts for the 
precession of the moving frame relative to the barycentric 
frame. We recall that the time coordinates are related 
by t = i+ c~ 2 A(t) + 0(c~ 4 ). The quantities A(t) and 
R a (t) that appear in the transformation are determined 
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by Eqs. ( 6.15 1 and ( 6.20 I , respectively. The inverse trans- 



formation is 



where 



£ ab {t) = N a c {t)N b d {t)S cd {t), 
B ab {t) = M a c (t)M b d (t)B cd (t), 

M ab {t) := 5 ab - \e abc R c (t) + 0^) 



(6.33) 



(6.34) 



is the inverse to M. ab - I n these equations we have t = 
t - c~ 2 A(t) + 0(c~ 4 ), and the quantities A(t) and R a (t) 
are determined by 



dA 
dt 



1 



-v + U ext , 



(6.35) 



^abc 



dR c 



-id [a U^-3v [a d b] U cxt . 



Expanding Eqs. (6.321 in powers of c 2 produces 
Sab = Sab + c- 2 [~Ad- t S ab + 2e cp{a RPS c b) ] + 0( C - 4 ) and 



B nb = B n 



0(c 2 ). In this we substitute Eqs. (6.25) 



(6.271, and we replace the hatted potentials by the orig- 



inal potentials using the correspondence of Eqs. (5.221 
and (5.231. After simplification we obtain 



S ab = -d ab U cxt + ^ {-d {ab) y cxt + 4v c (d ab U c c xt - d c{a Ulf) - id t{a U^ - 2(v 2 - U ext )d ab U eKt 

+3v c v {a d b)c U C xt + 2u( Q 3 b)t £/ 0Xt + 3d {a U cxt d b ) U cxt ] + 0(c" 4 ) (6.36) 



and 



B a b — ^ e pq(ad P b ^ (t^ex 



v"U cxt )+0(c- 2 ). (6.37) 



In these equations the external potentials are evaluated 
at x — z(t) after differentiation. Notice that unlike S ab , 
the barycentric tidal moment S ab does not involve the 
functions A and R a that must be obtained by integrating 
first-order differential equations; this was the reason for 



introducing the transformation of Eq. (6.32). Notice also 



that since B ab has been worked out to leading-order only 
its transformation to the barycentric frame is trivial. 



VII. TIDAL MOMENTS FOR A TWO-BODY 
SYSTEM 

The results obtained in the preceding section apply to 
any post-Newtonian environment described by the ex- 
ternal potentials C7 e xt! U^ xt , and ^ e xt- Given these po- 
tentials, the motion of the black hole i n the barycen- 
tric frame (t,x a ) is determined by Eq. (6.301, and the 



barycentric tidal moments S ab and B ab are obtained by 



evaluating Eqs. (6.36) and (6.371, respectively. The tidal 
moments perceived by the black hole are S ab and B ab , 
and these are calculated using the transformation of 



Eqs. (6.34l-(6.36l 



In this section we specialize our discussion to a spe- 
cific post-Newtonian environment that consists of a sin- 
gle external body (perhaps another black hole), so that 
the black hole is a member of a post-Newtonian two- 



body system. Adapting our notation to this specific sit- 
uation, we let our original black hole have a mass mi, 
position Zi(t), velocity i>i(i), acceleration ai(i), and so 
on. (These quantities were previously denoted m, z, v, 
and a, respectively. Below we will redefine m to be the 
system's total mass m\ + ni2, and v to be the system's 
relative velocity Vi —v%.) The second body also is mod- 
eled as a post-Newtonian monopole, and it has a mass 
m2 1 position Z2(t), velocity t>2(i), acceleration a-x (t), and 
so on. 

In Sec. VII A we list the potentials associated with the 
external body, evaluate their derivatives, and calculate 
the barycentric tidal moments. In Sec. VII B we simplify 
our expressions by writing them in terms of r :— z\ — 
22, the relative separation between the two bodies, from 
which the individual trajectories can be recovered. In 
Sec. VII C we restrict our attention to a binary system in 
circular motion, and in Sec. VII D we compute the tidal 
moments as viewed in the moving frame of the black hole. 
Finally, in Sec. VII E we compare our post-Newtonian 
answers to those obtained by Poisson [2] in the context 
of the small-hole approximation (see Sec. I A). 

Our calculations in this section rely on well-known re- 
sults from post-Newtonian theory. These can be found, 
for example, in Blanchet's review article [3"5] . 



A. Two-body potentials and tidal moments 

Assuming that the external body is a post-Newtonian 
monopole of mass m2, the external potentials can be 
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expressed in a form that is directly analogous to that 
of the black-hole potentials of Eqs. (14 lOh— (14 13h. We 



have C7 ext = Gm 2 /s, C/ c Q xt = Gm 2 v%/s, ip cxt = Gm 2 ^i 2 /s 
X oxt = Gm 2 s, and # ext = ipcx 1)2 



^d(X cyi t, where s now 



stands for \x — z 2 {t)\ and [i 2 := 



1, ,2 



Gmi/|zi - z 2 | 



These potentials are easily differentiated, and after 
evaluation at x = Z\(t) we obtain 



U ex t 

dabU ex t 
dtaUext 
dbUcxt 

dtbU^xt 



Gm 2 



r 

Gm 2 



Gm 2 



(3n a n b - S ab ), 



Gm2 (o x \ 

Gm 2 v 2 



Gm 2 v 2 



(3n b n c - 8 bc 



G 2 m 1 m 2 a 



n n b 



(7.1) 
(7.2) 
(7.3) 
(7.4) 
(7.5) 
(7.6) 



Gm 2 v% v c 
3 — {3n b n c - bc )v 2 , 



(7.7) 



dab^c 



Gm-, 



2vi 



Gmi 



(2>n a n b - 5 ab ) 



Gm 2 



2r 3 



3(n • v 2 ) 2 (5 ab - 5n a n b ) 

+ 12(n • v 2 )v 2 ( a nb) - 2v 2a v 2b 
G 2 m\m 2 



2r 4 



(Sab - n a n b ). 



(7.8) 



We have introduced the new quantities 

r:=zi-z 2 , r := \z\ - z 2 \, n:=r/r. (7.9) 



To arrive at Eqs. (7.7) and (7.8) we used the equations 
of motion a 2 = Gmin/r 2 + 0(c~ 2 ) to replace the ac- 
celeration vector of the second body by its Newtonian 
expression. 



Making the substitutions into Eqs. (6.361 and (6.371 

gives 



£ab 



3Gto9 



3Gm2 



l {ab) 



([2^-4^ 



172 ) + 2l>2 - -(TO ■ V 2 ) 



2 V ' 2 r 
- [3(to ■ i7i) - 2(n • 17 2 )] n {a v lb) + [4 (to • « x ) - 2 (to • v 2 )] n {a v 2b) 

+v 1{a v lb) - 2v 1{a v 2b) + v 2(a v 2b) j + 0(c" 4 ) 



2 5 Gmi Gttt-2 



l (af>) 



(7.10) 



and 



B. Generic orbital motion 



^ = -^J^ [» x («i - ^)] (a n 6) + 0(c- 2 ), (7.11) 

where := n a n b — ^S ab . The quantities A(t) 

and R a (t) that appear in the transformation from the 
barycentric frame to the black-hole frame are determined 
by the equations 



• 1 , Gm 2 9 . 
A = -v 2 + —^ + 0(c- 2 ), 



R = G ^nx{Av 2 -'iv l ) + 0{c- 2 ) ] 



(7.12) 



these are obtained from Eqs. (6.36) 



It is straightforward to generalize Eqs. ( 7.10 )-( 7.12) 



from a two-body system to an A^-body system by simply 
writing the external potentials as a sum of single-body 
terms. The generalized expressions can then be com- 
pared with the corresponding results of Damour, Sof- 
fel, and Xu — see, in particular, Eqs. (4.29)^(4.31) of 
Ref. jTTj . As was already stated in Sec. I B, we find that 
our expressions agree with theirs. 



To simplify the foregoing results we incorporate the 
fact that the motion of each body in a post-Newtonian 
two-body system can be related to the motion of the 
relative orbit, which is described by the separation vec- 
tor r = z\ — z 2 and the relative velocity vector v — 
v\ — v 2 . The post-Newtonian dynamics implies that 
if the system's barycenter is placed at the origin of 
the coordinate system, then Z\ = (m 2 /m)r + 0(c~ 2 ) 
and z 2 = —{mi/rri)r + 0(c~ 2 ), where m := mi + m 2 
is the total mass of the two-body system. As a con- 
sequence we also have V\ = (m 2 /m)v + O(o~ 2 ) and 
v 2 = —(m 1 /m)v + 0(c~ 2 ), and we make these substi- 



tutions in Eqs. ( 7. 10 )~( 7. 12 1 



In addition we incorporate the fact that the post- 
Newtonian motion of a two-body system takes place in 
a fixed orbital plane. We take this plane to be the x-y 
plane, and we use polar coordinates r and <p to describe 
the orbital motion. We write r = (rcos <j>, r sin </>, 0), and 
we resolve all vectors in the basis to = (cos 0, sin 0, 0), 
4> = (— sin 0, cos <j>, 0), and I = (0,0,1) associated with 
the polar coordinates; the vector I is normal to the plane, 
and it is aligned with the system's total angular momen- 
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turn. We have 



v = rn+ (r 



(7.13) 



and we also make these substitutions in Eqs. ( 7.10| |- 
(p2l. 



After simplification our results become 

£ab '- 



3Gm 2 3Gm 2 



3m 2 



2m 2 



2W>) 5 



5Gmi + 6Gm 2 



2r 



Sab 



(2mi + m 2 )m 2 ./ , 
^2 ?W)n (a </> 6) 

K^) 2 0<a O >) +0(c~ 4 ), 

6Gm 2 , , , _ 2s 
§— l( a n b ) + 0{c ), 



[r 2 + (rc, 



2m 2 

Gm 2 4mi + 3m 2 
2r 2 m 



Gm 2 



(r</>) l a + 0(c- z ) 



We recall that m = mi + m 2 is the total mass of the 
system, r := \z\ — z 2 | is the inter-body distance, <fr is 
the angular position of the relative orbit in the orbital 
plane, n is a unit vector that points from body 2 to 
body 1, <^> is a unit vector that points in the direction 
of increasing tfi, and finally, I is the unit normal to the 
orbital plane. We use the notation 7i( a b) = n a n b — |<5 ao 
and <j)( ab ) — <j) a <j>b- \&ab- 



C. Circular motion 

To specialize to circular orbits we set r = and r<p = v 
in the preceding results. The post-Newtonian equations 
of motion imply that 



cu :— — 



Gm 



l- l -{Z-r,){v/cY + 0{c-*) . (7.1M 



where m = mi + m 2 is total mass and r\ :— mim 2 /m 2 
is a dimensionless reduced mass. As a consequence of 
Eq. \7.18 \ we find that v = s/G m /r + 0(c~ 2 ). Making 
these substitutions in Eqs. ( 7. 14 )— ( 7.17 1 gives 



£ a b — — 



3Gm 2 



mi + 2m 2 2 

1 - -^2m— {v/c) 



l(ab) 



B ab 

A 

R a 



+ (V/C) 2 ^(ab)) +0(C" 4 ), 

-3— vl( a n b ) + G{c ), 

(2m x +3m 2 )m 2 2 _ 2 

2^ V + ° {C } ' 

Gm 2 Ami + 3m 2 



2r 2 



u/ a + 0(c" 2 ). 



(7.19) 
(7.20) 
(7.21) 
(7.22) 



For circular orbits we also have cj> = tot, with u given by 



Eq. (7.181 



To list the components of the tidal moments it is useful 
to define 



So 
£ic 
£is 

&2c 
£2s 



= ^11 + ^22), 

= £l3, 

= ^23 ! 

= ~ ^ 22 )' 



(7.23) 



and 





(7.14) 


Bo 




B22), 




(7.15) 


Bi c 
Bi s 


:= B13, 
:= B23, 






(7.16) 


B 2c 


■= \{B U - 


B22), 


- 2 ). 


(7.17) 


B 2s 


:= B12. 








With the vectorial 


basis n = 


(cos 0, sin <^,0), 



(7.24) 



(— sin 0, cos 0, 0), and I — (0,0,1) we find that the non- 
vanishing components of the tidal moments are 



S 

&2c 



Gm 2 
' 2r 3 
3Gm 2 
2r 3 



1 



1 - 



mi 
2m 
3mi + 4m 2 
2m 



{v/cf + 0(c- 4 ) 
{v/cY 



Bis 



+ 0(c" 4 ) 
3Gm 2 

+ 0(c" 4 ) 
3Gm 2 



cos 20, 

3mi + 4m 2 



2m 



(«/c) 5 



sin 20, 



., u cos + 0(c ), 
3Gm 2 



r 



: w sin + 0(c 



(7.25) 

(7.26) 

(7.27) 
(7.28) 
(7.29) 



The components £i Cl Ei s , Bo, i3 2c , and i3 2s all vanish 
for circular orbits. These results were already displayed 



in Sec. I A; in Eqs. ( 1.10 1— ( 1.14) we used the symbol b 



(instead of r) for the inter-body distance, w rc i (instead of 
v) for the relative orbital velocity, m (instead of mi) for 
the black-hole mass, and m! (instead of m 2 ) for the mass 
of the external body. 



D. Tidal moments in the black- hole frame 



The tidal moments of Eqs. ( |7.25| -( |7.29| ) refer to the 
barycentric frame. We may express them in the moving 
frame of the black hole by invoking the transformation of 



Eqs. (6.34 |-( 6.36 1. The transformation involves a switch 
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from global time t to local time t, and a rotation of the 
Cartesian axes mediated by the vector R(t). According 



to Eqs. (|7.21| the transformation of the time coordinate 
is given 



t = 



1 



(2mi + 3m2)m2 
2m? 



(v/c) 2 + 0(c~ A ) 



t. 



(7.30) 



And according to Eq. ( |7.22| we have c~ 2 i? a = -(flt)l a , 
with 



Gm (4m i + 3ra2)m2 
r 3 2m 2 



(v/c) 2 + 0(c- A ) (7.31) 



denoting the processional angular velocity of the moving 
frame relative to the barycentric frame. (This is the ro- 
tation of the coordinate axes, not the rotational motion 
of the black hole on its orbit.) The rotation takes place 
around the z-axis, and it is easy to show that it is effected 
by the transformation <\> — > <f> = <fi — ilt. 

Altogether we find that the tidal moments are given 

by 



Bo 
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Sin — — 
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3Gm 2 
2T 3- 

+ o( c - 4 ) 
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4:m 2 



2m 



B lc 



2r 3 

+ o( c - 4 ) 

3Gm 2 



1 - 



cos 20, 
3mi - 



4m; 



2m 



-(v/cY 



(v/cf 



sin 2(j>, 
wcos(/>+ Q(c~ 2 ), 



3Gms 



v sin 4> + 0(c 2 ) 



(7.32) 
(7.33) 

(7.34) 

(7.35) 
(7.36) 



in the moving frame. Tensorial express ions for £ a b and 
B a b can be obtained directly from Eqs. (7.191 and (7.20) 
by m aking t he sub stitut ion <fi — > </>. After involvement of 
Eqs. (17.181, (|7.30|, and (|7.3l| we find that <f> = Cot, with 



Gh 



1 



1 



(3 + r,)(«/c) 2 + 0(c- 4 ) 



(7.37) 



This is the angular frequency of the tidal moments 
as measured in the moving frame of the black hole. 
The transformation from w to w involves a switch from 
barycenter time to local proper time, and a rotation of 
the local accelerated frame relative to the global iner- 
tial frame. Notice the c hange in si gn in front of r\ := 
m\m 2 /m 2 between Eqs. (7.181 and (7.371. 



E. Comparison with Schwarzschild tidal fields 

The tidal moments of a black hole of small mass m\ 
moving in the gravitational field of another black hole 



of large mass m 2 can be obtained simply by evaluat- 
ing the components of the Riemann tensor for the large 
black hole; the Riemann tensor is evaluated in the mov- 
ing frame of the small black hole. The details of such 
a computation are presented in Poisson |14j . and in this 
subsection we compare our results. Poisson uses differ- 
ent definitions for the harmonic components of the tidal 
moments, and his results are presented in Schwarzschild 
coordinates. With the conventions adopted here, in rel- 
ativist's units, and in harmonic coordinates, Poisson's 
results are 



£o = 

&2c — 

&2s — 

Bic = 

B u = 



m 2 



2(r + m 2 ) 2 (r — 2m 2 ) ' 
3m 2 r — m 2 



2(r + m 2 ) 3 r — 2m 2 
3m 2 r — m 2 
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o 3/2 , 

6m 2 y/r — m 2 
(r + m 2 ) 3 r — 2m 2 

o 3/2 / 

6m 2 \Jr — m 2 
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cos 2</>, 
sin 2<j), 



COS ( 



sin ( 



where 



tot, with 



m 2 



(r + m 2 ) 3 



Here, r is the orbital radius of the small black hole (in 
harmonic coordinates), and t is proper time on the cir- 
cular orbit. 

It is easy to check that when (v/c) 2 :— m 2 /r -C 1, the 



Schwarzschild expressions reduce to Eqs. (7.32 )— ( 7.37l 
when mi <^ m 2 ; in this limit m ~ m 2 and 77 ~ 0. In 
their common domain of validity, our results agree with 
those of Poisson. 



VIII. GEOMETRY OF THE EVENT HORIZON 

In this section we present an application of the results 
obtained in Sec. VII. We examine, in a particular gauge, 
the intrinsic geometry of the event horizon of a tidally- 
deformed black hole. We emphasize that the discussion 
presented here is tied to a specific choice of gauge; a 
different slicing of the event horizon would produce a 
different intrinsic geometry. 

The harmonic coordinates (t, x a ) are singular on the 
black-hole horizon, and an examination of its geome- 
try requires a change of coordinates. For this pur- 
pose we return to the light-cone coordinates (v, p, 6 ) 
of Sec. Ill B. It is known [T5] that in the light-cone 
gauge, the coordinate description of the event horizon 
is p = 2Mi = 2Gmi/c 2 , the same as in the unperturbed 
Schwarzschild geometry. Equations (3.10 1— ( 3. 13 1 then 



imply that the induced metric on the event horizon is 
given by g AB = (2M 1 ) 2 fl AB + h AB + 0(Mi'R~ 2 C~ 1 ), 
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where £Iab is the metric on the unit two-sphere, and 

h AB = -~(2M 1 )*(£ AB + B AB ) (8.1) 

is the tidal perturbation. Here £ AB an d B\ R are the tidal 
potentials defined in Eqs. (3.16) and (3.181, respectively. 



As in Sec. Ill B, we (momentarily) refrain from displaying 
the factors of c as well as the overbar. 

To simplify the horizon metric we implement a gauge 
transformation generated by the vector field 



•2) 



where £\ and B A arc introduced in Eqs . (3.15) and (3.171, 
respectively. The transformation changes the metric per- 
turbation according to \iab — * h' AB — Hab — Da£,b — 
Db^Ai where Da is the covariant-derivative operator 
compatible with £Iab- Using the relations Da£ b — 
D R £! = ~ %Sl AB £* and D A B% 



J BC A — 2 C AB t^AB<^ ' cuiu ua^b 

we find that the new perturbation is given by 



D B B\ 



&AB> 



h' 



1 



AB 



(2M 1 ) 4 n AB £\ 



(8.3) 



where £ q is defined by Eq. (3.141 



Reinstating the factors of c and the overbar on the 
tidal moments (to emphasize that we are working in the 
black-hole's comoving frame), we find that the induced 
metric on the black-hole horizon is given by 



9ab = {2M 1 f 



1 — - ^ £ a b^ a ^ b 



AB 



o 



K 2 £ 
(8-4) 

in this choice of gauge. The line element on the horizon 
is 



ds 2 = (2M 1 f 



1 — - — - ^} £ a b£l a £l^ 



2c 2 



(d9 2 



<d<j> 2 ) 



(8.5) 



and £l a — (sin0cos0,sin0sin</>, cos0) is the unit radial 
vector. According to these equations, the area of each 
cross-section v = constant of the event horizon is given 
by A — 47r(2Mi) 2 , so that 2Mi is the radius averaged 
over each cross-section. 
The metric of Eqs. 



4 1 and (8.5 1 can be reproduced 



by embedding a closed two-surface described by 
r = 2M 1 [l + e(6,4>)] 



(8.6) 



in a flat, three-dimensional space charted by spherical 
coordinates (r, 0,</>). Working consistently to first order 
in e, we find that the metric on this two-surface is given 

„2 _ rtl.M2/i i Or\Ua2 



by ds 2 = (2M 1 ) 2 (l + 2e)(d9 2 
with Eq. (18. 51) when 



-sin Odtjr), and this agrees 



Ml 



£ a i,fl a D, b . 



•7) 



This equation evidently describes a quadrupole deforma- 
tion of a round two-sphere of radius 2M\ . To first order 
in e, the deformation produces no change in area. 

So far our considerations have been limited only by 
the restriction M\ -c TZ, which ensures that the tidal 
perturbation is small. As was discussed in Sec. I A, this 
restriction includes both the small-hole and weak-field 
approximations as limiting cases. If we now restrict our 
attention to the weak-field approximation and place the 
black hole within a post-Newtonian environment, then 
the tidal moments £ a & can be imported from Sec. VI D 
and substituted within Eq. (8.7). To illustrate this, we 



restrict our attention further to the situation examined 
in Sec. VII C, in which the black hole is a member of a 
two-body system in circular motion. The relevant tidal 
moments are listed in Eqs. ( 7.32 1— ( [7.34 1 , and when these 
are inserted within Eq. (8.7), we obtain 



M 1 2 M 2 
2b 3 



1 



Mi 
2M 



K cl /c) 2 + 0(c" 4 ) 



(1 - 3 cos 2 0) 



3MfM 2 
+ 2b 3 

+ 0(c" 4 ) 



3Mi + 4M 2 
2M 

9 cos 2ip. 



(frel/c) 5 



Here M\ := Gnii/c 2 is the black hole's gravitational ra- 
dius, and Mi := Gm-ijc 2 measures the mass of the com- 
panion body. We use the same notation as in Sec. I B: 
b is the separation between the two bodies (in harmonic 
coordinates), M := Mi + M 2 is a measure of the total 
mass within the system, and v rc i = cy 1 Mjb is the relative 
orbital velocity. The symbol tp stands for <f> — ujv, where 
v is the advanced-time coordinate on the event horizon 



and oj is the orbital frequency of Eq. (7.371 



Equation ( 8.8 1 implies that the event horizon is bulging 
along an axis directed toward the orbiting body. To see 
this clearly, we calculate from the metric the circumfer- 
ence of a line of longitude tjj = constant, and we obtain 



Ci = 



2tt(2Mi)|1 - 
3M 2 M 2 



46 3 



+ 0(c~ 4 ) 



1 



M(M 2 
4b 3 
3Mi 



Mi 
+ 2M 
4M 2 



2M 



cos 2tp 



( Vl . cl /c) 2 + 0(c~ 4 ) 
(v IC \/c) 2 

(8.9) 



This equation reveals that the circumference is largest 
(stretched) when ip — {0, 7r} and smallest (squeezed) 
when i/j = , ^r}. We also calculate the circumference 
of the equator (at = 5) and obtain 



Co 



2tt(2Mi) jl 
+ 0(c" 4 ) 



M\M 2 



1 



Mi 
2M 



(v IC \/cY 



(8.10) 



This equation also reveals a bulging of the horizon at the 
equator. 
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IX. TIDAL HEATING 

In this final section, we present another application 
of the results obtained in Sec. VII. We calculate the 
tidal heating of a black hole of mass mi placed in a 
post-Newtonian tidal environment created by an exter- 
nal body of mass 7712- For simplicity, we restrict our 
attention to circular motion. These results, unlike those 
presented in Sec. VIII, are gauge-invariant. In Sec. IX 
A we calculate the tidal heating of a nonrotating black 
hole, and in Sec. IX B we examine the case of a rotating 
black hole. The foundations for this calculation are given 
in Ref. [3D]. 



A. Nonrotating black hole 



in Sec. I B; in Eq. (1.171 we used the symbol b (instead 



of r) for the orbital separation, v re \ (instead of v) for the 
relative orbital velocity, m (instead of mi) for the black- 
hole mass, and m! (instead of 7712) for the mass of the 
external body. 



Equation (9.4 1 can be compared with the result ob- 
tained by Poisson [30 for a black hole of small mass m,\ 
moving in the field of another black hole of large mass 
m 2 . In geometrized units, and in harmonic coordinates, 
Poisson's result is 



mi 



32 /mi 

5 Vm.2 



'/m-2 
V r 



1 — 7712/V 



(1 + m 2 /r) 9 (l - 2m 2 /r) 5 



When m 2 /r = (v/c) 2 is small the relativistic factor be- 
comes 1 — 6(v/c) 2 + 0(c~ 4 ), and this expression agrees 
with Eq. (9.4) when mi -C m 2 . 



The rate at which a black hole of mass mi acquires 
mass by tidal heating is given by |30] 



Gmi 



16 (Gmi) e 



45 



,,15 



\B ab B ab +0(c- 4 ) 



(9.1) 



in which an overdot indicates differentiation with respect 
to t. This equation excludes contributions from octupole 
and higher-order tidal moments — see Ref. jTB] . It is easy 
to show, however, that these contributions occur at order 
c~ 4 (and smaller) relative to the dominant, quadrupolc 
term; they are therefore included in the neglected terms 
of Eq. dO. 



According to our results in Sec. VII D, the tidal mo- 
ments are given by 



S a b — — 



B ab = 



1 - 



3Gwi2 
r 3 

+ (v/c) 2 (f) {ab) ^ 
6Grri2 r 



mi + 2m,2 
2m 

+ 0(c- 4 ), 



(v/cf 



l (ab) 



(a n b) 



0(c 



(9.2) 
(9.3) 



in the black- hole frame. Here n a — (cos <f>, sin <f), 0) , 
a = (- sin 0, cos 0,0), and l a = (0,0,1), with <f> = 
Cot; the angular frequency to is displayed in Eq. (7.37|. 



The time derivatives of the basis vectors are given by 
n a = uxj) ai (f> a = —uin a , and l a — 0; this gives rise to 



l (ab) 



— 2u> n 



(afb) 



and 



(ab) = -2^n (a b) . 



Evaluating Eq. (9.1) fromEqs. (9.2 1 and (9.3 1 produces 
Gmi = 



32 m?mi/Gm\9 



5c 

5m 2 + Ylm\mi 



\ r J 
6m| 



1 



(v/c) 2 + 0(c- 4 ) 



. (9-4) 



where m = mi + m 2 is the total mass, r is the orbital 
separation (in harmonic coordinates), and v — U Gm/r 
is the relative orbital velocity. The rate at which the 
tidal coupling increases the black hole's angular momen- 
tum can next be obtained fr om the rigid-rotation rela- 
tion mic 2 = uiJ\. Equation (9.4 1 was already displayed 



B. Rotating black hole 

We next calculate the tidal heating of a rotating black 
hole, assuming that the tidal fields are not affected (at 
first post-Newtonian order) when the nonrotating black 
hole is replaced by a rapidly rotating hole. The rate at 
which the hole's angular momentum is increased by the 
tidal coupling is given by |30j 



GJi = - 



2 (GrV 5 



45 



„io 



-X^l + Sx^^i + c- 2 ^) 

3(4+17 X 2 )(£ 2 + c- 2 B 2 ) 
15 X 2 (E 3 + c- 2 B 3 ) + 0(c- 4 ) 



(9.5) 



where x := cJi/(Gm 2 ) is a dimcnsionlcss angular- 
momentum parameter that ranges between and 1, and 
E x := £ ab £ a \ E 2 := {£ ab s b ){E a c s c ), E 3 := (£ ab s a s b ) 2 , 
B, := B ab B ab , B 2 := (B ab s b )(B a c s c ), B 3 := (B ab s a S b ) 2 . 
Here, the unit vector s a points in the direction of the 
hole's spin angular-momentum vector, so that J\ = J\S. 
In this application, the spin and orbital angular momenta 
are aligned or antialigned, so that s a = el a with e = ±1. 
Evaluation of Eq. (9.5 1 produces 
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2 , mfm.2 (Gm) 7 



1 
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12 + 51x 2 m 2 
4 + 12x 2 m 



(v/cf 



+ 0(c- 4 ) 



(9.6) 



where 



mi + m 2 is the total mass, r the orbital 



separation (in harmonic coordinates), and v — U Gm/r 
is the orbital velocity. The rate at which the black-hole 
mass changes as a result of tidal heating can next be 
obtained from the rigid-rotation relation mic 2 = QJ\. 
With Eq. (|7.37| we get 



Gmi = -^r^x(l + 3x 2 )' 

DC 



777 
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12 X 2 
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12 X 2 
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(9.7) 



where the parameter e = ±1 was previously defined by 
the relation s° = el a . Thus, the black hole loses mass 
when the orbital motion proceeds in the same direction 
as the spinning motion (e = 1), and it gains mass when 
the orbital motion proceeds in the opposite direction (e = 
— 1). In each case the orbital motion is slower than the 
spinning motion, and the black hole always loses angular 
momentum. 



Equation (9.7) can be compared with the result ob- 



tained by Poisson j3U] for a rotating black hole of small 
mass mi moving in the field of a nonrotating black hole 
of large mass m 2 . In geometrized units, and in harmonic 
coordinates, Poisson 's result is 



mi 



-X(l + 3x 



15/2 



Y— 

5 ^ v A ' V?772 ) \ r 

15 2 

(1 -m 2 /r)(l- 4+1 | x2 m 2 /r) 
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When m2/r = (w/c) 2 is small the relativistic factor be- 
comes 



18 + 69 X 2 , , . a , _ 4 x 

1 TTT2F (w/c) +0(c } ' 



and this expression agrees with Eq. (9.7 1 when mi <C 7772. 



The same conclusion holds when the small black hole 
moves in the field of a large rotating black hole. In this 
situation, the error term in the previous expression is of 
order c -3 instead of order c~ 4 . 
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